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ABSTRACT 

The INT Photometric Ha Survey of the Northern Galactic Plane (IPHAS) is a 1800 
square degrees CCD survey of the northern Milky Way spanning the latitude range 
—5° < b < +5° and reaching down to r' ~ 20 (lOer). It may increase the number 
of known northern emission line sources by an order of magnitude. Representative 
observations and an assessment of point-source data from IPHAS, now underway, are 
presented. The data obtained are Wide Field Camera images in Ha alpha narrow-band, 
and Sloan r' and i' broad-band, filters. We simulate IPHAS (r' — Ha,r' — i') point- 
source colours using a spectrophotomctric library of stellar spectra and available filter 
transmission profiles: this gives expected colours for (i) solar-metallicity stars, without 
Ha emission, and (ii) emission line stars. Comparisons with Aquila field observations 
show that simulated normal star colours reproduce the data well for spectral types 
earlier than M. Spectroscopic follow-up of a Cepheus field confirms that sources lying 
above the main stellar locus in the (r' — Ha, r' — i') plane are emission line objects, 
with very few failures. Examples of Ha deficit objects - a white dwarf and a carbon 
star - are shown to be readily distinguished by their IPHAS colours. The role IPHAS 
can play in studies of nebulae is discussed briefly, and illustrated by a continuum- 
subtracted mosaic image of the SNR, Shajn 147. The final catalogue of IPHAS point 
sources will contain photometry on about 80 million objects, (abridged for astro-ph)) 
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1 INTRODUCTION 

The astronomical significance of Ha spectral line emission is 
that it both traces diffuse ionized nebulae and is commonly 
prominent in the spectra of pre- and post-main-sequence 
stars and binaries. Since these are objects in relatively short- 
lived phases of evolution, they are a minority in a mature 
galaxy like our own. Their scarcity has in turn acted as a 
brake on our understanding of these crucial evolutionary 
stages that in youth help shape the growth of planetary 
systems, and in old age determine stellar end states and 
the recycling of energy and chemically-enriched matter back 
into the galactic environment. The major groups of emission 
line stars include all evolved massive stars (supergiants, lu- 
minous blue variables, Wolf-Rayet stars, various types of 
Be star), post-AGB stars, pre-main-sequence stars at all 
masses, active stars and interacting binaries. This last group 
most likely harbours SN la progenitors within it, in guises 
that are still subject to considerable debate (Hillebrandt & 
Niemeyer 2000; Uenishi, Nomoto & Hachisu 2003). 

Existing catalogues of emission line objects contain any- 
thing from a few, to a few hundred, sources. Within the least 
populous object classes (e.g. the luminous blue variables and 
supersoft X-ray binaries, with just a few of each known in 
the Galaxy) there can be a confusing melee of 'special cases' 
that inhibit confident identification of essential and general 
behaviours. In effect, stellar evolutionary studies have been 
bedevilled by small number statistics and a lack of good de- 
mographics. The remedy for this problem is to exploit the 
technical developments of recent years that have boosted 
both the efficiency with which large scale astronomical sur- 
veys can be performed, and the quality that can be achieved. 
In particular, large area CCD-mosaic detectors offering good 
spatial resolution have now completely supplanted the pho- 
tographic techniques of the last century. In this paper, we 
describe the Isaac Newton Telescope (INT) Photometric Ha 
Survey of the Northern Galactic Plane (IPHAS), a pro- 
gramme that began taking data with the INT Wide Field 
Camera in the second half of 2003. 

The goal of IPHAS is to survey the entire northern 
Galactic Plane in the latitude range —5° < b < +5° - a 
sky area of 1800 sq.deg. The choice of latitude range was 
tensioned between the rising total telescope time require- 
ment and the expected fall off in discoveries to be made with 
increasing latitude (see below). The 10-degree wide strip re- 
quires in the region of 22 weeks clear time, and the hope is to 
complete the observations before the end of 2006. The data 
obtained will be mined both for spatially-resolved nebulae 
and for unresolved emission line stars. For point sources, the 
magnitude range will be 13 < r' < 20. Here we will focus 
on presenting the basic features of the survey, together with 
the extraction of point source data and the analysis of pho- 
tometric colour information. The different technical issues 
relating to the identification and measurement of resolved 
Ha-emitting nebulae will be presented in a later paper. For 
now, we just point to the opportunity that IPHAS presents 
both for making new discoveries and for high quality Ha 
emission mapping on large angular scales. 

To place this new northern hemisphere survey in con- 
text, it is appropriate to review the scale and character of the 
emission line star population that previous Galactic Ho sur- 
veys have revealed. Kohoutek & Wehmeyer (1999, hereafter 



KW99) have added their own discoveries within the latitude 
range -10° < b < +10° (1979 objects; data obtained in the 
years 1964-1970) to those of a wide range of independent 
searches: these go back as far as the original work of Mer- 
rill & Burwell that resulted in the Mount Wilson Catalogue 
(MWC, see Merrill & Burwell 1933). The total number of 
KW99 objects is 4174. In many cases the source observa- 
tions are spectra obtained using objective prism facilities. 
For sources in the northern hemisphere, this compilation 
supercedes that due to Wackerling (1970). Three-quarters 
of the stars listed in KW99 are assigned a photovisual mag- 
nitude, m pv < 13, and it is surmised that this is roughly the 
catalogue's completeness limit. They also note that over 80 
percent of all the objects they list in the —10° < b < +10° 
band fall within the narrower —5° < b < +5° band. At 
the fainter magnitudes we are exploring, we might expect 
this concentration toward the Galactic Equator to become 
even more pronounced. Naive extrapolation of the bright- 
end (m pv < 13) magnitude distribution of the KW99 emis- 
sion line stars to span 13 < m pv < 20 would suggest that 
our survey should uncover 8000-10000 new objects. 

This is probably an underestimate for a number of rea- 
sons. First, we can check this extrapolation of KW99 against 
the same quantity derived from the Stephenson & Sand- 
uleak (1971, hereafter SS71) southern Galactic Plane survey. 
The SS71 completeness limit is shallower at m pg ~ 11. We 
find, even on excluding the Galactic Bulge region located 
exclusively in the southern sky, that the prediction rises to 
~40000. This dramatic difference can have a number of ori- 
gins - beginning with simple differences in the Galactic stel- 
lar populations accessible from the northern and southern 
hemispheres, and ending with issues of experimental tech- 
nique. Nevertheless a parallel between the KW99 and SS71 
catalogues is that the bright magnitudes sampled strongly 
favour early-type, intrinsically luminous stars (such objects 
account for three-quarters of the KW99 catalogue). On going 
to much fainter magnitudes the sampled emission line star 
population is likely to broaden in character as intrinsically 
fainter object types (e.g. young and active stars, interacting 
binaries) become included. 

An immediate precursor to IPHAS and, indeed, a 
prompt for the need for a northern survey, is the 
AAO/UKST narrow-band Ha Survey of the Southern 
Galactic Plane and Magellanic Clouds. This was the last 
photographic sky survey carried out on the UK Schmidt 
Telescope (UKST). It was completed in 2003 and is now 
available as digital survey data derived from SuperCOSMOS 
scans of the original survey films (the SHS database, located 
at http://www-wfau.roe.ac.uk/sss/halpha/). A description 
of this survey is presented by Parker et al (2005) : important 
points to note are its high spatial resolution (~ 1 arcsec) and 
its areal completeness - the entire southern Galactic Plane 
was imaged within the latitude range —10° < b < +10°. 
Each of the 233 Galactic Plane fields observed had an effec- 
tive dimension projected on the sky of 4° x 4°. The south- 
ern survey has provided the source material for a variety 
of continuing research projects (see e.g. Morgan, Parker & 
Russeil 2001; Parker & Morgan 2003; Drew et al 2004). For 
the detection of point sources, IPHAS betters both the SHS 
sensitivity and spatial resolution, and offers the advantage 
of CCD dynamic range and linearity. The sensitivity of the 
two surveys to spatially-resolved Ha emission is comparable. 
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Table 1. Co-ordinates of centres, estimates of maximum Galactic reddenings, observation dates and seeing for the IPHAS fields discussed 
in this paper. Note that the offset partner field (nnnno to field nnnn) is offset by 5 arcmin N and 5 arcmin E. The reddening estimates 
are derived from the Schlegel et al (1998) reddening maps. 



IPHAS Field Centre Coordinates maximum observation mean 



field 


RA 


Dec 




b 


reddening 


date 


seeing 




number 


(2000) 


(2000) 


(°) 


(°) 


[E(B - V)] 


(dd/mm/yyyy) 


(arcsec) 




2540 


05 33 49 


+25 15 


181.73 


-4.18 


1.1 


05/11/2003 


0.8 


Section 2, 5 


4090 


18 32 03 


+00 41 


31.33 


+4.62 


2.3 


09/06/2004 


1.1 


Section 3 


4095 


18 33 33 


+01 47 


32.48 


+4.79 


1.6 


09/06/2004 


1.0 


Section 3, 4 


4199 


18 47 27 


+01 58 


34.23 


+1.78 


3.0 


12/06/2004 


0.9 


Section 3, 4 


6985 


22 14 56 


+61 11 


105.18 


+3.83 


2.3 


05/11/2003 


0.8 


Section 6 


6993 


22 15 56 


+61 44 


105.59 


+4.21 


2.5 


05/11/2003 


0.8 


Section 6 


7012 


22 18 36 


+61 22 


105.65 


+3.73 


1.8 


03/11/2003 


1.2 


Section 6 


7019 


22 19 39 


+61 55 


106.06 


+4.12 


2.4 


03/11/2003 


1.2 


Section 6 



We begin the description of IPHAS in the next section 
with a presentation of the niters used, and our observing 
and data reduction techniques. Following this, in section 3, 
we discuss the use of the Ha,r',i' filter photometry in the 
diagnostic (r' — Ha) versus (r' — i') diagrams that can be 
constructed from the survey data for point sources. Specif- 
ically, we introduce simulated colour-colour tracks for both 
normal solar-metallicity stars occupying the main stellar lo- 
cus and for emission line objects. We then provide examples 
of (r — Ha,r' — i') diagrams in three contrasting northern 
Galactic plane locations (sections 4 to 6). The fields dis- 
cussed are identified in Table 1. In section 4, we illustrate 
the application of the simulated tracks for normal stars with 
reference to fields in Aquila; in section 5 we perform a con- 
sistency check of IPHAS photometry of a Taurus field ob- 
tained on a photometric night; and in section 6 we present 
some follow-up spectroscopy relating to a field in Cepheus 
that illustrates the high success rate achieved in the con- 
firmation of candidate emission line objects. In Section 7, 
we outline the application of IPHAS to imaging spatially- 
resolved nebulae, and illustrate this with the beautiful ex- 
ample of the supernova remnant, S 147. The paper ends with 
a summarising discussion (Section 8). 



2 SURVEY OBSERVATIONS AND DATA 
EXTRACTION 

2.1 IPHAS observations 

The Wide Field Camera, mounted on the 2.5-metre Isaac 
Newton Telescope, is an imager comprising 4 AR-coated, 
thinned 4K x 2K EEV CCDs arranged in an L shape, cap- 
turing data from an on-sky area of approximately 0.3 of 
a square degree. With a pixel dimension of 13.5 /im, cor- 
responding on-sky to 0.333x0.333 arcsec 2 , the instrument 
is appropriately configured to fully exploit the high quality 
sub-arcsecond seeing frequently encountered at the Roque 
de los Muchachos Observatory in La Palma. Adequately 
sampled ~1 arcsec resolution is particularly useful given that 
lower-reddening Galactic Plane star fields, observed down to 
~20th magnitude, are at times very crowded. 

Not accounting for the geometric consequences of the 
L-shaped detector arrangement, the total number of point- 
ings required to span the 10 x 180 deg 2 survey area would be 
6000. On accounting for the detector outline and requiring 



a little overlap between pointings, we have chosen to fix the 
number of field centres at a total of 7635. Furthermore, each 
pointing is paired with a second pointing at an offset of 5 
arcmin W and 5 arcmin S, such that the number of quality- 
controlled sets of exposures expected to be compiled into 
the final survey database is 15270. Stars falling into a gap 
between the mosaiced CCDs in one exposure are captured 
in a partner exposure set. Nevertheless the great majority 
of Galactic Plane sources will be imaged at least twice. The 
pairing and offsetting, together with the chosen tessellation, 
comes very close to complete coverage of the northern Plane 
(> 99 %). 

Since Ha falls in the red part of the spectrum, IPHAS 
was conceived of as a large-scale programme that could 
readily make use of less heavily subscribed bright and grey 
nights. To ensure this scheduling flexibility, whilst obtaining 
the associated continuum-band observations required for es- 
tablishing unambiguous Ha excesses, it was decided to re- 
strict our broadband choices to red or longer wavelengths. 
This also has the effect of increasing the penetration of the 
survey for a given exposure time since these longer wave- 
lengths are also less subject to Galactic dust obscuration 
than UBV bands. This stands in contrast to the somewhat 
bluer emphases of the older Ha catalogues (e.g. SS71, where 
the objective prism data spanned 3300 < A(A)< 6800). The 
particular choice we made was to obtain the Ha exposures 
alongside Sloan r' and i' filter observations. 

All three filter profiles are plotted as Fig. 1. The Sloan 
niters have been preferred over Harris alternatives because of 
their squarer transmission profiles. The Sloan r' filter is the 
most blue-sensitive of the three (central wavelength 6240 A), 
with the Ha filter positioned toward the red end of its band- 
pass (central wavelength 6568 A). With a full- width at half- 
maximum (FWHM) transmission of 95 A, the Ha filter is 
more than broad enough to capture all likely Doppler shifts 
due to Galactic motions of up to a few hundred km s _1 or 
~10 A, as well as blueshifts of up to an additional ~10 A 
due to the converging beam of the INT/WFC. The cen- 
tral wavelength of the Sloan i' filter is 7743 A. We have 
added two broad band filters to this survey in order to give 
a continuum-dominated colour with which the (r' — Ha) ex- 
cess measurement can be compared. It has been shown in the 
past (see e.g. Robertson & Jordan 1989) that this is impor- 
tant for distinguishing between a genuine Ha emission excess 
and a molecular band dominated late type stellar spectrum 
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Figure 1. The transmission profiles of the Ha, Sloan r 1 and 
i' filters used in all IPHAS observations. The r' and i' fil- 
ter profiles are plotted as solid lines, whilst the Ha profile is 
shown dashed. The dotted line is the mean WFC CCD re- 
sponse at the cooled working temperature. These profiles are 
available in numerical form from the INT/WFC web page, 
http://www.ing.iac.es~quality/filter/filt4.html 

- in such cases large (r' — Ha) 'colour' will correlate with 
relatively extreme (r' — i'). In truth, the diagnostic value of 
this strategy is wider than this, as shall become apparent in 
Section 3. 

The exposure times in the three filters were set at 
120 sec (Ha) and 10 sec (r' and i') for the first season's 
observing in 2003. Evaluation of these data, once extracted, 
led us to increase the r' band exposure to 30 sees, from the 
start of the 2004 observing season, to compensate better for 
their typically higher moonlit background. This adjustment 
also acknowledges the pivotal role the r' band exposures 
must play in the survey's exploitation - it is important that 
errors in this band, appearing in both the Ha excess and the 
broadband colour measurement, are minimised. 

For the purpose of photometric calibration, each night's 
observations includes standard fields, obtained in twilight 
and at intervals of approximately 2 hours through the 
night. The standards are chosen from a list including the 
Landolt equatorial fields (Landolt 1992), Sloan (Smith et 
al 2002) and Stetson standards (at the Canadian Astron- 
omy Data Centre, http://cadcwww.dao.nrc.ca/standards/). 
Nightly observations are also acquired of spectrophotomet- 
ric standards with a view to assisting the final calibration 
of the narrow-band Ha data. A programme of supporting 
spectrophotometric observations is planned with a view to 
placing the Ha calibration on the desired firm footing in the 
longer term. 

2.2 Data processing 

Processing of IPHAS INT WFC data generally follows 
the pipeline procedure devised by Irwin and Lewis (2001) 
for dealing with optical mosaic camera data. The two- 
dimensional instrumental signature removal includes provi- 
sion for: non-linearity correction at the detector level; bias 
and overscan correction prior to trimming to the active de- 



tector areas; flatfielding; and fringe removal in the i! pass- 
band. 

Flatfielding in all bands is accomplished by stacking 
suitable twilight Hatfield exposures taken over the course of 
each typically one-week observing run to create master cal- 
ibration flats. These have been found to be stable on this 
timescale provided no filter changes, or other instrumental 
setup changes, occur in the middle of the run. The gain 
differences between each detector in each passband are re- 
moved by normalising a robust measure of the average sky 
level for each detector to a common system (in this case the 
sky level on CCD no. 1). The flatfielded i' data, even for the 
short exposures (10 s) used here, show measurable fringing. 
At the same time, the data taken for the IPHAS project are 
not themselves suitable (short exposures in crowded Galac- 
tic Plane regions) to construct good quality fringe maps for 
correcting this problem. Since the fringing in the i' band 
is relatively stable with time, we make use of a library of 
i' band fringe maps taken from other observing runs using 
the INT WFC. These have been found to reduce the level of 
fringing to an acceptable level when used with the defringing 
algorithm in the pipeline. 

Each master flat, in conjunction with a previously de- 
fined bad column list, is also used to construct confidence 
maps for each passband. These are used during the cata- 
logue generation to flag less reliable pixels in each image by 
providing a measure of the inverse variance weight for each 
pixel e.g. bad pixels have zero weight, heavily vignetted re- 
gions have low weight, poor DQE pixels have lower weight, 
and so on. These confidence measures are used directly to 
weight the image detection part of the catalogue generation 
algorithm and help avoid generating excessive numbers of 
spurious images around defects and other excessively noisy 
regions. 

Catalogue generation follows the precepts outlined by 
Irwin (1985, 1997) and includes the facility to: automati- 
cally track any background variations on scales of typically 
20-30 arcsec; detect and deblend images or groups of im- 
ages; and parameterise the detected images to give various 
(soft-edged) aperture fluxes, position and shape measures. 
The generated catalogues start with an approximate World 
Coordinate System (WCS) defined by the known telescope 
and camera properties (eg. WCS distortion model) and are 
then progressively refined using all-sky astrometric cata- 
logues (eg. USNO, APM, 2MASS) to give internal precision 
generally better than 0.1 arcsec and global external preci- 
sion of 0.25 arcsec with respect to USNO and APM, and 
0.1 arcsec with respect to 2MASS. These latter numbers are 
solely dependent on the accuracy of the astrometric cata- 
logues used in the refinement. 

All catalogues for all CCDs for each pointing are then 
processed using the image shape parameters for morpholog- 
ical classification in the main categories: stellar; non-stellar; 
noise-like. A sampled curve-of-growth for each detected ob- 
ject is derived from a series of aperture flux measures as a 
function of radius. The classification is then based on com- 
paring the curve-of-growth of the flux for each detected ob- 
ject with the well-defined curve-of-growth for the general 
stellar locus. This latter is a direct measure of the integral 
of the point spread function (PSF) out to various radii and 
is independent of magnitude, if the data are properly lin- 
earised, and if saturated images are excluded. The average 
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Figure 2. Magnitude and colour errors in IPHAS data as a func- 
tion of magnitude. The quantity shown, for each filter in the top 
3 panels, is the rms deviation per magnitude bin, between the 
measured magnitudes for each of IPHAS fields 2540 and 2540o. 
The lower two panels show the colour rms deviations. The bin 
size is 0.5 magnitudes for all but the i' data, where a binning 
of 0.25 magnitudes was more convenient. See further discussion 
of fields 2540/2540o in section 5. These data were obtained on a 
photometric night in November 2003. The r' exposure times were 
10 sec in this case. 

stellar locus on each detector is clearly defined and is used 
as the basis for a null hypothesis stellar test for use in classi- 
fication. The curve-of-growth for stellar images is also used 
to automatically estimate frame-based aperture corrections 
for conversion to total flux. 1 

Any photometric standards observed during the run 
(mainly Landolt 1992 and spectrophotmetric standards) are 
automatically located in a standards database and used to 
estimate the zero-point in each passband for every pointing 
containing any of these standards. The trend in the derived 
zero-points is then used to assign a photometric quality in- 
dex for each night and also as a first pass estimate for the 
magnitude calibration for all the observations. The Ha fil- 
ter is treated as equivalent to a standard Johnson- Cousins 
R-band filter to obtain a Vega-like magnitude which is used 
as an initial calibration (to be refined later, as mentioned 
above) . 

Various quality control plots are generated by the 
pipeline and these are used to monitor characteristics, such 



We note that in regions of intense nebular emission with in- 
creasingly short spatial scale variations of the " background" , au- 
tomatic detection, parameterisation and classification of objects 
becomes progressively more unreliable. In such regions continuum 
subtraction via difference imaging will yield better results. 



as:- the seeing; the average stellar image ellipticity (to mea- 
sure trailing); the sky brightness and sky noise; the size of 
aperture correction for use with the "optimal" aperture flux 
estimates (here "optimal" refers to the well-known property 
that soft-edged apertures of roughly the average seeing ra- 
dius provide close to profile fit accuracy eg. Naylor 1998). 
The "optimal" catalogue fluxes for the r' , i' and Ha filters 
for each field are then combined to produce a single matched 
merged catalogue from which diagnostic colour-magnitude 
diagrams and two-colour diagrams may be produced. These 
merged catalogues - the fundamental IPHAS product - con- 
tain flux, classification and match position error for each 
object in each passband. 

The tAU-registered naming conven- 

tion for all point sources derived from these catalogues is 
IPHAS JHHMMSS.ss+DDMMSS.s - thereby encoding the 
2000 object co-ordinates into the name. 

To give an impression of the internal magnitude errors 
in the catalogued magnitudes and derived colours we plot, in 
Fig. 2, the rms deviation between the magnitudes measured 
in each filter, and the associated colours, for point sources 
common to two overlapping exposure sets (fields 2540 and 
2540o, discussed again in Section 5). These were obtained 
on a photometric night in November 2003 as the moon was 
setting. Calculated empirically as \/ < (m2540 — WI25400) 2 > 
or its colour equivalent, over a range in mean magnitude 
Am = 0.5 or Am = 0.25, the error is corrected back to a 
representative single field measurement error by dividing by 
V2. The bright-end errors in the magnitudes themselves, in 
plots such as these, are typically dominated by calibration 
offsets of a few hundredths that will be removed when a 
final uniform survey calibration is devised. In the case of 
fields 2540 and 2540o the offsets were all small (less than 
~ 0.01). From 2004 on, when the r' exposures were increased 
to 30 sec, the faint-end r' errors drop to around 60% of 
those for 2003 (for the same sky conditions). This carries 
through to the colour errors falling to 80% or less of their 
2003 levels. Altogether, this significantly raises the fraction 
of catalogued objects that will meet the quality target of 
Ar' < 0.1 for r' < 20. 

To date roughly 3 Tbytes of raw data from the first 
two seasons of IPHAS observing have been processed this 
way. This corresponds to well over 100,000 4kx2k CCD 
images and over 40 million objects have been catalogued. 
All of this processed data is also available at the individ- 
ual frame and catalogue level via a PostgreSQL database 
interface which allows users to: postage stamp browse for 
candidate verification; construct image catalogue overlays, 
including on-the-fly matching with other catalogues such 
as the 2MASS point source catalogue; perform on-demand 
continuum image subtraction and mosaicing of larger ar- 
eas; access all the quality control information; and more 
(see Irwin et al 2005). The database interface is available 
on the Cambridge Astronomical Survey Unit (CASU) web- 
site, at http://apm2.ast.cam.ac.uk/cgi-bin/wfs/dqc.cgi. Co- 
ordinates of the centres of the observed IPHAS fields are 
obtainable there. 2 



2 These are identified via object names taking the form int- 
phasjmnn*, where nnnn is a 4-digit number up to 7635, and * is 
the wild card for further characters identifying exposure type. 
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3 SIMULATION OF THE IPHAS 
COLOUR-COLOUR PLANE 

The three bandpasses of the survey provide the basis for 
the construction of a number of magnitude-colour diagrams 
and a colour-colour diagram to describe any chosen region in 
the northern Galactic Plane. Using just the two r' , i', broad 
bandpasses, one may derive colour-magnitude diagrams that 
can in principle reveal different sequences at different red- 
denings that may be present in the field under investigation. 

Full exploitation of IPHAS hinges on the colour-colour 
plane involving all three bands. The combination of magni- 
tudes we use is (r 1 — i') as abscissa and (r' — Ha) as ordi- 
nate, so that objects with Ha band excesses appear higher 
within the diagram, while intrinsically redder or more highly 
reddened objects are over to the right. The most straightfor- 
ward use that can be made of such diagrams is to pick out for 
spectroscopic follow-up those objects whose (r' — Ha) colour 
places them clearly above the main locus of non-emission line 
objects. Additional information contained within the colour- 
colour diagrams can lead to identification of more subtle 
candidate emission line stars and also to a characterisation 
of the stellar populations distributed along the line of sight. 
In this second sense, IPHAS can also be seen as providing a 
far-red map of stellar populations in the northern Galactic 
Plane. 

Fig. 3 is a composite (r 1 — Ha,r' — i') plot derived 
from data obtained in three paired IPHAS fields (fields 
4090/4090o, 4095/4095O and 4199/4I99o: see Table 1). In 
each case, catalogues of sources classified by the CASU 
pipeline as either 'definitely stellar' or 'probably stellar' were 
extracted from within a 30 x 30 arcmin 2 box spanning most 
of the overlap region between the two pointings. For each ex- 
tracted object, the datum plotted is the mean of the colours 
derived independently from each of the two exposures mak- 
ing up the field pair. The data shown are limited to the 
magnitude range 13 < r' < 20, where the error in either 
colour is kept to less than ~0.05 magnitudes. These are rep- 
resentative of the better data in the IPHAS database in that 
they were obtained on photometric nights in June 2004 at 
times of ~1 arcsec seeing and low sky background. 

All three fields are located in the Aquila Rift region, 
and sample sightlines that pass through the outer parts of 
the molecular cloud. Dame & Thaddeus (1985) noted that 
this is a nearby (~200 pc) and not particularly opaque cloud 
system, presenting around 2 magnitudes of visual extinction 
only. This is a modest addition to the reddening through the 
remaining Galaxy beyond - the reddening data of Schlegel, 
Finkbeiner & Davis (1998) indicate maximum visual extinc- 
tions, Av ranging from ~ 5, in 4095, up to ~ 10, in 4199. 
The nearby rift cloud is responsible for the lightly-populated 
gap, seen in Fig. 3, between the upper sequence and the 
lower, but much more densely populated strip. The exis- 
tence of this separation allows a clear demonstration of how 
well theoretically-synthesised tracks compare with and make 
sense of the photometry. 

To achieve an understanding of the behaviours seen in 
the colour-colour domain, we have constructed two types 
of synthetic tracks: the first type concerns the properties of 
normal stars without Ha emission, while the second explores 
the impact of adding narrow Ha emission to generic stel- 
lar spectral energy distributions (SEDs). We present these 
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Figure 3. A plot of (r' — Ha) versus (r' — i') for 13,818 point 
sources extracted from the overlap area of the paired exposures 
for IPHAS fields 4199 (black points), 4090 (red points) and 4095 
(blue points), in Aquila. The sources lie in the magnitude range 
13 < r' < 20. The colours for each object have been formed 
independently from each set of Ha, r' and i' exposures, and then 
averaged. The median error in either colour is ~ 0.03 at r' close to 
19, and only occasionally exceeds 0.06 at the r' = 20 limit of the 
plot. Note the very well defined upper edge to the upper sequence 
— this is where unreddened main sequence stars are located. The 
few stars above it are candidate emission line objects. 

tracks below, using the Aquila fields to illustrate the former 
in section 4. 



3.1 The IPHAS colours of normal stars 

For simulating the (r 1 — Ha) and (r' — i') colours of nor- 
mal stars, we have used the library of stellar spectral en- 
ergy distributions (SEDs) due to Pickles (1998, hereafter 
P98). At a final binning of 5 A the spectra in this li- 
brary are well enough sampled that we may use them to 
compute narrow-band Ha relative magnitudes with confi- 
dence, alongside the analogous broadband r' and i' quan- 
tities. The required numerical filter transmission profiles, 
shown in Fig. 1, are available via the ING WFC web pages 
(http:/ /www.ing. iac.es/Astronomy/instruments/wfc/), as 
is a mean Wide Field Camera CCD response curve. To en- 
sure compliance with the Vega-based zero magnitude scale, 
we have defined synthetic colour as follows: 

(r - O = -2.5 log( — - - AA ) + 2.5 log( ^-^- )(1) 

where T' r and T[ are the r' and i' numerical transmission 
profiles, after multiplying by the mean WFC CCD response 
curve, and rebinning to match the P98 spectral library sam- 
pling. The SED for Vega, F\ t v, is the appropriately resam- 
pled version of that due to Hayes (1985). The (r' — Ha) 
colour is evaluated in the same way, after substituting the 
Ha numerical profile in place of the i' profile. Since Vega 
is an A0V star, its SED at Ha incorporates a strong ab- 
sorption line feature. Currently, because the CASU pipeline 
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uses broad-band standard fields to calibrate measured source 
magnitudes, there is an offset in (r' — Ha) colour between 
the catalogue data and our simulations. The dominant spec- 
tral type in the standard fields will be appreciably later 
than Vega's AO, with the consequence that the standard- 
star SEDs will both be redder and less eroded by Ha line 
absorption. On this basis one would expect, and we do find, 
that zero (r' — Ha) for unreddened main sequence stars cor- 
responds to (r'—i 1 ) ~ 0.3 (late F), rather than to (r' — i') = 
(Vega, AOV), in plots of IPHAS data obtained in photomet- 
ric conditions. Hence, on comparing simulated tracks with 
observation it is necessary to correct for this. The best way 
to do this is to assume that only a shift in (r 1 — Ha) is 
required, given that the calibration of the broadband-only 
(r — i') colour should be secure enough (see Section 5 and 
the discussion of Fig. 9). The shift that needs to be applied 
to theoretical (r — Ha) values to match them to observa- 
tion is then always downwards, varying in amount between 
about -0.10 and -0.25. 

We have simulated tracks for main sequence stars (lumi- 
nosity class V) , giant stars (class III) and supergiants (class 
I) using, for simplicity only, solar metallicity P98 spectra. 
Each sequence has been calculated for a range of redden- 
ings using an R = 3.1 optical-IR extinction law in the form 
given by Howarth (1983). The colours derived are given for 
Eb-v = 0, 1, 2, 3 and 4 in Tables 2 and 3. They are speci- 
fied at this level of detail because there is no single redden- 
ing vector that translates all of a sequence directly onto its 
reddened counterpart. The unreddened dwarf, giant and su- 
pergiant sequences are compared in Fig. 4, where it can be 
seen that there is a gradual decrease in track gradient with 
increasing luminosity class. Nevertheless between mid-F and 
late-K spectral types there is minimal distinction between 
the luminosity classes - this is the regime where Ha absorp- 
tion is weak and there is not yet any marked development 
of the molecular band structure that typifies M-type stars. 

The effect of interstellar extinction on the sequences is 
illustrated in Fig. 5. With increasing reddening, the range of 
(r' — Ha) colour spanned by each luminosity-class sequence 
diminishes. This shrinkage in dynamic range is due to effec- 
tive wavelength of the r' bandpass lengthening, and moving 
closer to the Ha bandpass as reddening becomes ever more 
extreme. It may be seen in Fig. 5 that the supergiant se- 
quence essentially reddens along itself, while the main and 
giant sequences shift across in a manner that sweeps out 
area in the colour-colour plane. Also shown in this figure 
is the reddening locus for AO dwarfs - for non-degenerate 
stars and binaries this line amounts to an important bound- 
ary in that none should drop below it. In practice, stars will 
be carried into the forbidden domain by a range of types of 
observational error. Degenerate dwarfs and related objects 
with stronger Ha absorption than early A stars would also 
fall below this line. 

3.2 The IPHAS colours of emission line stars 

To assess the impact on (r' — Ha) and (r' — i') colours of 
increasingly strong Ha emission, we have represented un- 
derlying stellar SEDs using simple power laws or, at later 
spectral types, blackbodies. At the acceptable price of some 
approximation, such as ignoring stronger spectral features 
like the Paschen limit, this approach allows us to explore a 



broad range of SEDs flexibly and to quantify thresholds for 
the straightforward detection of Ha emission. We present 
results for 4 simple SEDs: three power laws of the form 
Fx oc A - * 3 with 13 set equal to 4 (Rayleigh- Jeans case, rele- 
vant to the hottest O stars), 3 (appropriate to ~A0 stars) 
and 2.3 (the optically thick accretion disk case); a Planck 
function at a temperature of 5900 K that is a good match 
to the G2V SED in the P98 library. 

Ha emission, where present, takes on a wide range of 
profiles in stellar spectra, and can be practically any width 
- with FWHM anywhere in the range from a few 10s of km 
s _1 up to 1000s. But, for now, we treat the simple limiting 
case of Ha emission that is well- contained within the width 
of the INT/WFC narrow-band Ha filter. The particular re- 
alisation used is of a rectangular profile of breadth 25 A 
centred at 6570 A (note that the effective Ha filter band- 
pass is blueshifted for objects observed off-axis, which is why 
the central wavelength of the WFC filter is accordingly spec- 
ified as 6568 A). This yields results negligibly different from 
using a Gaussian profile. 

The synthesis of colours consists of the following steps: 
an underlying stellar SED is chosen; a rectangular Ha emis- 
sion profile of the desired equivalent width (EW) is super- 
imposed; the resultant artificial spectrum is then reddened 
as required using the reddening law specified in section 3.1; 
finally, the reddened SED is multiplied by the product of 
the survey filter profiles and WFC response and integrated 
to form colours as in equation 1. We have not attempted 
to apply this procedure to very late type SEDs dominated 
by molecular bands - in these stars, neither can the SED 
be easily parameterised, nor is an objective definition of Ha 
EW straightforward. 

On the basis of this procedure we have synthesised 
colours for the same set of reddenings (E(B — V) = to 
4 in steps of 1) for each of the 4 adopted stellar SEDs. Our 
results are presented in Fig. 6 and in Table 4. We find that 
there is a practical degeneracy between reddening and un- 
derlying SED such that a very-nearly unique locus is traced 
at each adopted Ha emission EW. This means that, in prin- 
ciple, a given location in the colour-colour plane, above the 
main stellar locus, is associated with a particular Ha EW. 

A further property of the SED-specified tracks is that 
at EW up to ~100 A, the trend with increasing EW is 
nearly vertical. But as Ha EW becomes very large, the 
tracks bend toward smaller (r' — i') as the Ha emission be- 
comes a more significant contributor to the r' flux. Indeed 
beyond an EW of 1000 A as the switch from a 'stellar' to a 
'nebular' spectrum with little discernable continuum takes 
place, the bending becomes very extreme. In reality (r' — i') 
in the nebular case will also depend somewhat on the rel- 
ative strength of line emission in the i' band - left out of 
consideration here. The limiting value of (r' — Ha) in the 
absence of any continuum is ~ 3.24 for our synthetic system 
referred to the Hayes (1985) SED for Vega (see Table 4). 
At the present time, without a properly defined zero point 
to the Ha filter magnitudes, this translates to an effective 
observed upper limit on (r' — Ha) of around 3.1 - any value 
appreciably above this signals a problem with the individual 
object's photometry. 

Finally an important feature to note in the trend in 
IPHAS colours with respect to both emission EW and red- 
dening is that the threshold for the detection of Ha emission 
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Figure 5. The effect of interstellar extinction, calculated according to an R = 3.1 Galactic law, on synthetic stellar tracks in the 
(r' — Ha,r' — i') plane. As in Fig. 4, the main sequences are drawn in black, the giant sequences in red, and supergiants in green. The 
three sets shown apply to Eb—v = 0.0, 2.0 and 4.0, as labelled. The dashed line shown is the reddening locus for A0V stars. This defines 
a notional minimum line for all non-degenerate stars. It is named the "early-A reddening line" for reasons that become clear in Section 
5. 
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Table 2. Synthetic tracks in the (r' — Ha, r' — i') plane, for main-sequence dwarfs and giants, calculated for a range of reddenings. 



Spectral Reddening:- 



Type 


EjB—Y — 


0.0 




= 1.0 


Eb — v - 


= 2.0 


^B — V 


= 3.0 


^B — V 


= 4.0 




W - i') 


(r' - Ha 


) (r' - i' 


) (r' - Ha) 


(r' - i') 


(r' - Ha) 


(r' - i' 


) (r' - Ha) 


(r' - i' 


) ( r ' - - 


05V 


-0.176 


0.067 


0.493 


0.274 


1.146 


0.438 


1.784 


0.560 


2.408 


0.643 


09V 


-0.151 


0.072 


0.519 


0.278 


1.173 


0.441 


1.811 


0.563 


2.436 


0.645 


BOV 


-0.143 


0.057 


0.529 


0.264 


1.184 


0.427 


1.824 


0.549 


2.451 


0.632 


B1V 


-0.120 


0.069 


0.552 


0.273 


1.209 


0.435 


1.850 


0.554 


2.478 


0.635 


B3V 


-0.049 


0.056 


0.623 


0.258 


1.279 


0.418 


1.919 


0.535 


2.546 


0.614 


B8V 


-0.022 


0.059 


0.649 


0.259 


1.304 


0.416 


1.944 


0.532 


2.572 


0.609 


B9V 


-0.016 


0.045 


0.652 


0.243 


1.305 


0.399 


1.943 


0.514 


2.569 


0.590 


AOV 


0.029 


0.001 


0.699 


0.199 


1.352 


0.355 


1.991 


0.468 


2.616 


0.544 


A2V 


0.046 


0.023 


0.716 


0.220 


1.370 


0.374 


2.009 


0.486 


2.634 


0.560 


A3V 


0.068 


0.025 


0.737 


0.220 


1.391 


0.372 


2.029 


0.483 


2.654 


0.555 


A5V 


0.107 


0.035 


0.776 


0.228 


1.428 


0.379 


2.065 


0.488 


2.689 


0.558 


A7V 


0.138 


0.070 


0.805 


0.261 


1.457 


0.408 


2.094 


0.515 


2.717 


0.584 


FOV 


0.212 


0.114 


0.876 


0.299 


1.525 


0.441 


2.159 


0.542 


2.780 


0.606 


F2V 


0.257 


0.156 


0.920 


0.338 


1.567 


0.477 


2.201 


0.576 


2.822 


0.638 


F5V 


0.263 


0.185 


0.926 


0.364 


1.574 


0.502 


2.208 


0.600 


2.830 


0.661 


F8V 


0.325 


0.218 


0.985 


0.393 


1.631 


0.526 


2.262 


0.620 


2.882 


0.677 


GOV 


0.368 


0.204 


1.028 


0.378 


1.673 


0.510 


2.303 


0.602 


2.921 


0.658 


G2V 


0.384 


0.227 


1.043 


0.400 


1.687 


0.530 


2.317 


0.621 


2.935 


0.676 


G5V 


0.394 


0.239 


1.052 


0.411 


1.695 


0.540 


2.323 


0.630 


2.941 


0.684 


G8V 


0.420 


0.261 


1.077 


0.429 


1.718 


0.556 


2.345 


0.643 


2.961 


0.695 


KOV 


0.445 


0.278 


1.098 


0.443 


1.736 


0.567 


2.360 


0.651 


2.973 


0.700 


K2V 


0.474 


0.298 


1.127 


0.461 


1.765 


0.582 


2.390 


0.665 


3.004 


0.713 


K3V 


0.548 


0.303 


1.199 


0.463 


1.836 


0.581 


2.460 


0.660 


3.073 


0.705 


K4V 


0.613 


0.347 


1.262 


0.502 


1.895 


0.615 


2.515 


0.691 


3.126 


0.732 


K5V 


0.650 


0.390 


1.301 


0.544 


1.937 


0.658 


2.561 


0.735 


3.175 


0.778 


K7V 


0.800 


0.490 


1.450 


0.636 


2.085 


0.743 


2.709 


0.812 


3.322 


0.849 


MOV 


0.903 


0.499 


1.553 


0.641 


2.188 


0.743 


2.811 


0.807 


3.424 


0.839 


M1V 


1.063 


0.547 


1.719 


0.688 


2.360 


0.789 


2.987 


0.853 


3.604 


0.885 


M2V 


1.144 


0.624 


1.795 


0.756 


2.431 


0.849 


3.055 


0.907 


3.669 


0.934 


M3V 


1.521 


0.820 


2.174 


0.941 


2.811 


1.023 


3.436 


1.070 


4.052 


1.088 


M4V 


1.829 


0.889 


2.470 


0.995 


3.096 


1.061 


3.710 


1.095 


4.315 


1.099 


M6V 


2.514 


1.060 


3.165 


1.143 


3.801 


1.191 


4.425 


1.207 


5.041 


1.198 


08III 


-0.180 


0.075 


0.489 


0.282 


1.142 


0.446 


1.780 


0.568 


2.404 


0.651 


B1-2III 


-0.108 


0.072 


0.565 


0.276 


1.222 


0.438 


1.863 


0.558 


2.491 


0.639 


B3III 


-0.042 


0.089 


0.630 


0.291 


1.285 


0.449 


1.925 


0.566 


2.551 


0.644 


B5III 


-0.019 


0.083 


0.653 


0.283 


1.308 


0.441 


1.947 


0.557 


2.573 


0.634 


B9III 


-0.022 


0.031 


0.649 


0.231 


1.303 


0.388 


1.942 


0.504 


2.568 


0.582 


AOIII 


0.013 


0.022 


0.682 


0.220 


1.335 


0.375 


1.973 


0.489 


2.598 


0.565 


A3III 


0.041 


0.057 


0.708 


0.252 


1.358 


0.404 


1.994 


0.515 


2.617 


0.587 


A5III 


0.107 


0.089 


0.772 


0.280 


1.420 


0.429 


2.053 


0.537 


2.674 


0.606 


A7III 


0.130 


0.096 


0.797 


0.285 


1.447 


0.432 


2.083 


0.539 


2.706 


0.607 


FOIII 


0.177 


0.103 


0.844 


0.290 


1.497 


0.435 


2.134 


0.539 


2.759 


0.606 


F2III 


0.230 


0.188 


0.890 


0.368 


1.534 


0.506 


2.164 


0.603 


2.782 


0.664 


F5III 


0.269 


0.177 


0.930 


0.356 


1.575 


0.492 


2.206 


0.588 


2.825 


0.648 


GOIII 


0.378 


0.236 


1.035 


0.408 


1.677 


0.537 


2.306 


0.627 


2.923 


0.681 


G5III 


0.449 


0.279 


1.102 


0.444 


1.739 


0.567 


2.364 


0.651 


2.977 


0.700 


G8III 


0.474 


0.276 


1.128 


0.440 


1.767 


0.563 


2.392 


0.647 


3.007 


0.696 


KOIII 


0.505 


0.284 


1.157 


0.445 


1.793 


0.564 


2.417 


0.646 


3.029 


0.692 


K1III 


0.529 


0.300 


1.182 


0.461 


1.820 


0.580 


2.445 


0.660 


3.059 


0.706 


K2III 


0.574 


0.325 


1.225 


0.481 


1.860 


0.596 


2.483 


0.673 


3.095 


0.716 


K3III 


0.582 


0.320 


1.230 


0.475 


1.863 


0.589 


2.483 


0.665 


3.092 


0.706 


K4III 


0.670 


0.351 


1.320 


0.502 


1.955 


0.611 


2.578 


0.683 


3.191 


0.721 


K5III 


0.808 


0.381 


1.456 


0.525 


2.088 


0.629 


2.708 


0.695 


3.318 


0.728 


MOIII 


0.935 


0.473 


1.581 


0.611 


2.212 


0.708 


2.830 


0.768 


3.438 


0.796 


M1III 


0.984 


0.461 


1.633 


0.598 


2.266 


0.695 


2.887 


0.755 


3.496 


0.782 


M2III 


1.040 


0.512 


1.679 


0.641 


2.304 


0.730 


2.916 


0.782 


3.519 


0.803 


M3III 


1.311 


0.564 


1.957 


0.687 


2.586 


0.770 


3.203 


0.817 


3.809 


0.832 


M4III 


1.652 


0.644 


2.296 


0.753 


2.924 


0.822 


3.539 


0.857 


4.144 


0.860 


M5III 


1.941 


0.725 


2.581 


0.821 


3.205 


0.879 


3.816 


0.904 


4.417 


0.899 
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Table 3. Synthetic tracks in the (r' — Ha, r' — i') plane for supergiant stars calculated for a range of reddenings. 



Spectral Reddening:- 



Type 


Eb — V — 

(V - i') 


0.0 

{r 1 - Ha) 


Eb — V — 

rv - %') 


1.0 

(r' - Ha) 


Eb — V — 

(r' - %') 


2.0 
(r' - Ha) 


Eb — V — 

rv - %') 


3.0 
(r' - Ha) 


Eb—v — 
(r' - %') 


4.0 

(V - Ha) 


BOI 


-0.078 


0.109 


0.591 


0.311 


1.243 


0.470 


1.880 


0.588 


2.504 


0.666 


B3I 


-0.063 


0.126 


0.603 


0.325 


1.253 


0.481 


1.888 


0.596 


2.510 


0.672 


B5I 


0.029 


0.147 


0.693 


0.341 


1.340 


0.491 


1.973 


0.601 


2.594 


0.672 


B8I 


0.049 


0.157 


0.717 


0.351 


1.369 


0.502 


2.007 


0.612 


2.631 


0.683 


AOI 


0.060 


0.135 


0.726 


0.327 


1.377 


0.477 


2.012 


0.586 


2.636 


0.657 


A2I 


0.071 


0.154 


0.734 


0.345 


1.382 


0.493 


2.015 


0.601 


2.636 


0.671 


FOI 


0.191 


0.131 


0.856 


0.318 


1.505 


0.462 


2.140 


0.565 


2.762 


0.631 


F5I 


0.219 


0.158 


0.882 


0.340 


1.529 


0.480 


2.162 


0.580 


2.783 


0.643 


F8I 


0.248 


0.183 


0.908 


0.363 


1.553 


0.500 


2.184 


0.597 


2.803 


0.657 


GOI 


0.317 


0.229 


0.971 


0.401 


1.609 


0.531 


2.234 


0.622 


2.847 


0.677 


G2I 


0.357 


0.241 


1.010 


0.412 


1.648 


0.540 


2.272 


0.629 


2.885 


0.682 


G5I 


0.389 


0.282 


1.040 


0.449 


1.676 


0.573 


2.298 


0.658 


2.909 


0.708 


G8I 


0.471 


0.312 


1.109 


0.468 


1.731 


0.584 


2.340 


0.661 


2.940 


0.704 


K2I 


0.611 


0.329 


1.251 


0.481 


1.876 


0.591 


2.488 


0.663 


3.089 


0.701 


K3I 


0.680 
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Figure 6. The effect of adding in increasing Ho emission to a range of simplified stellar SEDs. The black lines, representing an unreddened 
main sequence and the early-A reddening trend, are the expected bounds to the main stellar locus. The dashed lines are lines of constant 
Ha emission equivalent width, with the label at the left specifying the equivalent width in Angstroms. The vertical coloured lines are, in 
effect, Ha emission curves of growth for particular choices of underlying SED and reddening: the darker blue lines are for the Rayleigh 
Jeans case (F\ oc A~ 4 ) at different reddenings; light blue and green lines are for power law indices of -3 and -2.3 respectively, while the 
red lines are the results for a 5900 K blackbody (a G2V star, roughly). The lowest point on each curve corresponds to a narrow net Ha 
absorption of EW, 10 A. The reddenings, as E(B — V), are specified across the top and refer to each set of coloured lines. 



is lowest for bluer and/or less reddened objects (see Fig. 6). 
This implies, for instance, that IPHAS will pick out faint, 
nearby accreting objects very well indeed down to just a few 
Angstroms EW. Conversely, in the worst case of a densely 
populated main stellar locus spanning a wide range of red- 
denings, the EW threshold on the straightfoward detection 
of classical T Tau stars at E(B - V) ~ 2 [Av ~ 6) is 
around 30 A. Not infrequently, however, at larger (r' — i') 
(j> 2) the colour-colour plane below the unreddened main se- 



quence may be sparsely populated - the few objects located 
here could be the result of a combination of anomalous red- 
dening and line emission. Indeed it is generally the case that 
objects, checked as having reliable photometry, lying outside 
the bounds of the densely populated main stellar locus for 
their field, have a relatively high probability of being inter- 
esting in one way or another. 
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Table 4. Synthetic tracks in the (r' — Ha, r' — i') plane, for increasing Ha emission equivalent width. Data are provided for four 
underlying simplified stellar SEDs (3 power laws and a blackbody) and for a range of reddenings. The usual sign convention for Ha EW 
(left hand column) is reversed in that positive values refer to net emission. 



Ha Reddening:- 



EW (A) 
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E B -v = 1.0 
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1.043 
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1.125 
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0.991 




1.348 


1.619 




1.460 


2.237 


1.536 


200 


-V.ooo 


1.1 1 ( 


0.317 


1.365 


0.953 




1.512 


1.577 




1.620 


2.192 


1.693 


300 


-0.386 


1.441 
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2.232 
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3.242 
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" 3 (~A0 SED) 
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2 3 (optically thick disk accretion) 
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4. 5900 K blackbody (~G2 V SED) 
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Figure 7. The (r'—Hct, r' —i') plane derived from observations of 
IPHAS fiels 4095/4095o, compared with selected synthetic tracks. 
In both panels, the uppermost solid line traces the unreddcncd 
MS track, while the dashed line is the early-A reddening line. See 
text for specification of other tracks. In the upper panel, green 
points indicate r' < 16, while blue indicates r' > 16. Below, 
yellow corresponds to 19.5 < r' < 20, green to 19 < r' < 19.5 
and blue to 18 < r' < 19. 

4 SIMULATED AND OBSERVED IPHAS 
COLOUR-COLOUR DIAGRAMS 
COMPARED - FIELDS IN AQUILA 

With the assistance of the synthetic tracks derived in the 
preceding section it is possible to begin to make sense of the 
morphologies appearing in IPHAS (r' — Ha,r' — i') diagrams 
- with a view to their fuller exploitation. For this purpose 
we have selected some of the highest quality IPHAS obser- 
vations obtained from pointings in the Aquila Rift region, 
allowing us to exploit its distinctive and easily interpreted 
colour-colour domain morphologies. 

We begin with field 4095 in Aquila, included in Fig. 3 
as the blue data points. This field is roughly centred on 
I = 32.5°, b — +4.8°, sampling a region close to the survey's 
Galactic latitude upper limit. The reddening data of Schlegel 
et al (1998) indicates that Eb-v typically does not exceed 
1.6 in this direction. This makes it the least obscured of the 
three fields and, correspondingly, the field presenting the 
highest apparent density of stars (7097 out of the 13818 in 
Fig. 3). In Fig. 7 the data on extracted point sources are 
compared with selected synthetic tracks that have all been 



shifted downwards in (r' — Ha) by 0.17 to best match them 
to the data. 

The upper panel in the Fig. 7 shows the brighter end of 
the magnitude range (13 < r' < 18) that includes a modest 
number of very nearly unreddened M dwarfs and a much 
larger number of mainly giant stars. Indeed the M giants 
form a particularly tight sequence at (r' — i') 2.0. This 
suggests that most of the Galactic reddening along this sight 
line accumulates nearby because, if it were not, we would 
expect to see a more smeared giant distribution. The sim- 
ulated giant tracks for Eb-v = 1.4 and 1.6 are compared 
with this very well-defined feature. For (r' — i')<2 many of 
the brighter objects will be giants at a plausible reddening; 
but at (r' — i') > 2, the synthesised tracks for M2-5 III stars 
fall too low by ~ 0.05 in (r' — Ha). A similar problem affects 
comparisons between synthesised and observed tracks for M 
dwarfs also (see below). At r' < 18, only one object falls 
significantly below the early-A reddening line - it is likely to 
be a white dwarf or related object. 

The lower panel in Fig. 7 presents the faint end of the r' 
magnitude range, with the synthetic main sequence and gi- 
ant tracks, reddened to E(B — V) = 1.6, superimposed. The 
main locus of observed objects is now a little more steeply 
angled, indicating that these fainter stars include a much 
increased component of main sequence objects. However, at 
r' ~ 20, stars later in spectral type than mid-K are only de- 
tectable at E(B — V) < 0.8, as evidenced by the scatter of 
points extending the main stellar locus up to (r 1 — Ha) ~ 0.9. 
The small proportion of the plotted objects falling below the 
early-A reddening line can be presumed consistent with ob- 
servational error. The 0.17 offset of the synthesised tracks 
was determined by optimising the positioning of both this 
notional line and the unreddened main sequence with re- 
spect to the data for r' < 18. These particular IPHAS ob- 
servations have captured objects out to the limits of the 
Galactic disc population, such as reddened mid-M giants at 
~ 10 kpc, located around 800 pc above the mid-plane at 
about the location of the far Sagittarius-Carina arm. 

The more highly reddened Aquila field 4199, centred at 
I — 34.3°, b = +1.8°, provides some degree of contrast with 
4095 and is illustrated in Fig. 8. The reddening here is more 
variable with position, as well as more extreme. This shows 
itself directly in the broader red giant locus. Superimposed 
on the colour-colour plot are synthesised giant-star tracks for 
Eb-v — 2.4, and 3.0. The latter value was selected because 
the maximum Galactic extinction for this field, derived from 
the Schlegel et al 1998 mapping data, is Eb-v — 3.0. There 
is a rough consistency here with the findings from field 4095, 
in that the M giant track synthesised for this extinction falls 
a bit below the observed thinning of putative M giants (as 
it did for field 4095). Down to r = 20, the colour-colour 
data suggest the presence of main sequence stars of K and 
earlier type at reddenings in the range 1.2 E(B — V) 
2.4: the main sequence tracks for these limits are drawn to 
illustrate this (Fig. 8). This field again bears the imprint 
of the Aquila Rift in the relative deficit of stars between 
the lightly populated zero-extinction main sequence and the 
dense locus of stars at Eb-v > 1-2. 

In field 4199 it is more apparent that the observed 
unreddened M dwarfs tend to maintain the locus gradient 
defined at earlier spectral types, rather than begin to turn 
over as the synthesised track indicates they should. This 
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Figure 8. The (r' — /fa, r' — i') plane derived from observations 
of IPHAS fields 4f99/4199o, compared with selected synthetic 
tracks. The solid black lines are main sequence tracks for E(B — 
V) = and 2.4, while the dashed line is the early-A reddening 
line. The two red tracks are giant tracks, for E(B — V) = 2.4 and 
3.0. The position of the probable C star discussed in the text is 
picked out by the arrow. 

is most likely another symptom of the problem behind the 
M-giant discrepancy. At the present time, the available con- 
version between Landolt (R — I) colours, appropriate to the 
standard star fields, and Sloan (r —i') are not properly de- 
fined for M dwarf colours (see Smith et al 2002). Similarly, 
the existing conversion used in the CASU pipeline is not 
validated for (r' — i') > 1.5. Clearly this will need to be 
corrected in the future. In the mean time, the comparison 
between observations of M stars and synthesised data will 
be increasingly qualitative as (r' — i') increases beyond the 
validation limit. 

In neither field 4095 nor field 4199 do supergiants stand 
out in any obvious morphological way. This is likely to 
be both a consequence of their relative rarity and of the 
way in which their locus shifts almost along itself with in- 
creasing reddening. In principle, extremely red, isolated ob- 
jects located below the red giant locus could be picked out 
as candidate reddened supergiants - or as potential exam- 
ples of other interesting object types. Indeed, we find that 
IPHAS J184644.25+015324.6, the one isolated object in this 
part of the 4199 colour-colour plane at (r' — i',r' — Ha) = 
(3.36 ±0.02, 0.499 ±0.03), cannot be a reddened supergiant. 
One reason is that the reddening (E(B — V) ;> 4) required 
to explain its position in these terms is excessive relative 
to the maximum expected for the field. Another is that 
the 2MASS point source within 0.2 arcsec of this object's 
position exhibits very bright JHK magnitudes with un- 
usual colours (K = 9.55 ± 0.02, (J - H) = 1.68 ± 0.03, 
(H - K) = 0.76 ± 0.03). At r' = 18.75 ± 0.01, there 
are no grounds for doubting the reliability of the IPHAS 
photometry and the reality of the source. The absence of 
any significant proper motion rules this object out as a 
nearby brown dwarf. This object is largely absent from pre- 
existing photographic surveys, except that there is a de- 
tection of it in the UKST Infrared (IVN) Survey reported 
on the SuperCOSMOS Sky Survey website (http://www- 
wfau.roe.ac.uk/sss/): it is reported there at an / magnitude 
of 18.383, around 3 magnitudes fainter than the IPHAS i! 



magnitude. The IVN plate was obtained in 1981. This vari- 
ability, the anomalously low (r' — Ha) colour, together with 
its NIR colours, point towards a carbon star at a reddening 
corresponding to E(B — V) ~ 1.4 (see e.g. Bessell & Brett 
1988). A faint-end absolute K magnitude for such a star, 
if on the AGB, would be -6.5 (Claussen et al 1987). This 
places it at ~15 kpc. 



5 A COMPARISON BETWEEN 

FLUX-CALIBRATED SPECTRA AND IPHAS 
PHOTOMETRY - A FIELD IN TAURUS 

Up to this point the interpretation of the (r 1 — Ha,r' — i') 
plane has been based on synthetic photometry derived from 
P98 library spectra. Early in February 2004, we obtained 
WHT/ISIS service spectra of 6 stars selected from IPHAS 
data on fields 2540 and 2540o in Taurus. These IPHAS im- 
ages, obtained on 5th November 2003, were among the first 
to be pipeline-processed and were picked for closer investi- 
gation as examples of apparently good quality data obtained 
in good seeing and photometric conditions. The aim of the 
follow-up service spectra was to obtain relative spectropho- 
tometry in order to ascertain optical SEDs, spectral types 
and reddenings for the sample stars as a retrospective check 
on the IPHAS colours, and the typical errors in them. This 
exercise gives an impression both of the current state of the 
photometric calibration of the data and of the quality of the 
synthetic colour comparisons derived from P98. 

Field 2540 is centred on RA 05 33 49 Dec ±25 15 00 
(2000) in Taurus, only a few degrees from the Galactic an- 
ticentre direction. The characteristics of this sky position 
are very different from those in Aquila: here, the maximum 
Galactic reddening is modest and more smoothly varying, 
ranging from E(B - V) ~ 0.8 in the NE of the ~ 30 x 30 
arcmin 2 extracted region up to ~ 1.1 in its SW. The other 
obvious difference, which stands out in the (r' — Ha, r — i') 
plot for 2540/2540o in Fig. 9, is the absence of any red gi- 
ants within the magnitude range shown (13 < r' < 20). This 
absence is not just a consequence of the imposed magnitude 
limits since a K/M giant at 10 kpc viewed through ~3 visual 
magnitudes of extinction should be detected at r' ~ 17: it 
must be a real absence. In the example of solar and lower 
metallicity isochrones presented by Bertelli et al (1994), a 
red giant branch is only well-developed from around 100 
million years of age onwards - suggesting that the stellar 
populations sampled in this part of the outer Galaxy are 
younger than this. 

8 stars were initially selected for ISIS spectroscopy from 
the colour-colour diagram for fields 2540/2540o on the ba- 
sis that they were not too faint (r' < 18) and lay on the 
outer boundary of the main stellar locus (dark blue aster- 
isks in Fig. 9). These criteria biased the selection in favour 
of evolved spectral types. In the event, 6 of the 8 stars 
(which we refer to as stars A-F) were observed during a 
service night of mediocre weather. Both the blue and red 
arms of ISIS were used, with the R600B and R316R grat- 
ings installed, delivering spectra spanning 3500-5000 A and 
6000-8700 A for each star. To further the aim of relative 
spectrophotometry, the slit width was set fairly wide at 1.8 
arcsec, while the slit orientation tracked the parallactic an- 
gle. The resolution of the spectra is ~ 3.6 A. An observation 
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Table 5. Observed and derived properties of the follow-up sample of stars in IPHAS fields 2540/2540o in Taurus. 





IPHAS name/position 


IPHAS photometry 


spectral 


Eb-v 


synthetic 


colours 


comment 




J[RA(2000)+Dec(2000)] 


r' 


r'-i' 


r' - Ha 


type 




r'-i' 


r' - Ha 




A 


J053430. 11+251400.9 


13.19 


0.33 


-0.07 


A2V 


0.40 


0.32 


0.09 




13 


J053432. 14+252231.0 


15.82 


0.54 


0.12 


A0V 


0.65 


0.47 


0.26 


Ha emission 


C 


J053425. 99+250843.1 


13.01 


0.61 


0.16 


G0III 


0.35 


0.61 


0.30 




D 


J053458. 93+252316.7 


17.22 


1.91 


0.78 


M4V 


0.05 


1.86 


0.92 




E 


J053311. 04+251444.5 


15.36 


1.15 


0.30 


G0III 


1.10 


1.11 


0.42 




F 


J053305.64+251837.6 


17.06 


0.86 


0.13 


A5 


~1.1 


0.84 


0.25 


noisy ISIS data 
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Figure 9. The (r' — Ha,r' — i') diagram for IPHAS fields 2540 and 2540o, sampling an area roughly equivalent to 30x30 sq.arcmin 
around RA 05 33 40, Dec +25 20 00 (J2000) in Taurus. The magnitude range shown is 13 < r' < 20. The dark blue asterisks mark 
the colours of the 6 objects for which WHT/ISIS service spectra were obtained (stars A-F, see also Fig. 10 and Table 5). The red 
asterisks mark the 'predicted' colours derived from the flux-calibrated spectra, after applying the same shift to all 6 objects in order 
to minimise the mean difference with respect to the original IPHAS colours. Like the computed colours for stars A— F, the synthesised 
early-A reddening line has been shifted by +0.03 and —0.14 in (r' — i') and (r' — Ha), respectively. 



of the white dwarf G191— B2B was also obtained to serve as 
a spectrophotometric flux standard. 

The data were extracted from the CCD frames and then 
wavelength- and flux-calibrated using routines from the soft- 
ware package, FIGARO. The extracted 1-D spectra were 
then imported to the software tool, DIPSO, in order to de- 
termine approximate spectral types and reddenings by com- 
paring them with P98 library spectra. In every case, the 
spectral type determination rested on matching absorption 
line characteristics. This matching was performed using the 
blue spectra for all but the M4V star (star D) - for this ob- 
ject the red spectrum was more appropriate. For each star, 
the P98 library spectrum of the appropriate spectral type 
was progressively reddened, using the mean Galactic extinc- 
tion law (Howarth 1983), to identify the best fitting colour 
excess. The observed spectra and the best fits to them de- 
rived in this manner are shown in Fig. 10. The positions, 
magnitudes and further data on the 6 stars appear in Ta- 
ble 5. For star F, the data are not of sufficient quality to pro- 
vide a reasonable fix on luminosity class: but we suspect that 
its Ha profile indicates a lower gravity than a main sequence 



A5 star. In order to: smooth errors due to irregularities in 
the spectrophotometric flux calibration; avoid the need to 
correct for telluric absorption; provide a good extrapolation 
of the observed spectral energy distributions (SEDs) to cover 
the full spectral range, the final step of deriving photometric 
colours used the closest matching, appropriately-reddened 
P98 library spectra in the r' and i' bands, rather than the 
calibrated observations. Only the Ha fluxes were computed 
by multiplying the Ha filter profile directly with the cali- 
brated ISIS spectra. 

Both the original IPHAS colours and the colours derived 
from the fits to the spectrophotometry are listed in Table 5. 
As there is not yet a uniform and fully-verified zero-point 
calibration for all IPHAS frames, we have to shift the SED- 
based colours on to the IPHAS colours. The shifts that min- 
imise the mean differences in each of (r 1 — Ha) and (r' —i') 
are —0.14 and +0.03, with final rms deviations between the 
6 pairs of observed and predicted colours of 0.015 and 0.026, 
respectively. The size and sense of shift in (r' — Ha) is as 
expected (see section 3). That the shift in (r' — i') is small, 
but apparently finite, indicates that the night the IPHAS 
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Figure 10. Fits of P98 library spectra (green) to flux-calibrated WHT/ISIS spectra (black) of selected stars in IPHAS field 2540/2540o 
in Taurus. The label in each box (letters A to F) refer to the positions marked in the colour-colour diagram for this field (Fig. 9). Next 
to each label is the spectral type and colour excess (E(B — V)) estimate for the star. In each plot, the flux scale is linear, with zero as 
the minimum plotted value. The ISIS spectra are uncorrected for telluric absorption. 



imaging was obtained was not perfectly photometric. On 
the basis of the errors estimated for the IPHAS photometry 
(see Fig. 2) , we would expect the rms deviations between the 
catalogued and observed colours to be ~0.01. The somewhat 
larger values of 0.015 and 0.026 obtained here turn out to 
be determined mainly by errors in the relative spectropho- 
tometry and its analysis: for example, the uncertainty in the 
E(B — V) estimates is typically 0.05 and translates into a 
(r' — i') error of ~0.03. In (r' — Ha), the discrepancies are 
smaller and mainly arise in the r' band integration. The cir- 
cle from photometry to spectroscopy, back to photometry, 
closes satisfactorily. 

At the end of this process, it can be seen in Fig. 9 that 
star A (A2V) falls a little below the early-A reddening line, 
rather than just above, while the reddening line itself lies 



~ 0.02 magnitudes above the bottom edge of the main stel- 
lar locus. A part of the reason for this may be illustrated 
in Fig. 11 where the P98 A2V spectrum is superimposed on 
the ISIS observation of star A: it is possible that the limited 
resolution of the P98 library spectra (R ~ 500 or AA ~ 13 A 
at Ha) leads to the Ha in-band fluxes of early A-type stars 
being overestimated, very slightly. Another factor will be 
linked to the question of the mean Ha absorption EW and 
its variance for early-A stars as a function of sub-type and 
metallicity. The Ha absorption EW for star A is 14.5±0.4 A, 
while that for the P98 A2V star is 10.8 ± 0.5. These num- 
bers for A2V may be compared with the P98 A0V and Hayes 
(1985) Vega Ha EWs, that are both close to 13.0 A. Finally 
we note the contrast between the statement by Jaschek & 
Jaschek (1987) that the Balmer lines are strongest at A2 
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Figure 11. Excerpts of the spectra obtained for stars A and B, 
focusing on the He* line profiles. The WHT/ISIS observations are 
plotted in black, while the P98 library spectra of the appropriate 
spectral types (respectively A2V and A0V) are shown superim- 
posed in blue. The excess emission equivalent width in star B is 
12 A. In both panels the positioning and FWHM of the IPHAS 
Ha filter is indicated by the horizontal bar. 



(see their Table 10.1) and the maximum in P98 for near- 
solar metallicity dwarfs at AO. The underlying practical dif- 
ficulty here is the measurement and calibration challenge of 
the very broad Hcv absorption wings in A-type spectra. As 
the IPHAS survey completes and a uniform photometric cal- 
ibration is constructed, it will then be appropriate to sort 
the issue out and improve the absolute registration of this 
lower boundary. For the present, the A0V spectrum in the 
P98 library defines the shape of the early-A reddening line 
well enough to allow it to be used in a relative manner. 

The Ha profile for star B is also presented in Fig. 11. 
This object shows a distinct central reversal in Ha which 
allows it to be described as a weak emission line object. 
The excess (emission) equivalent width with respect to the 
library A0V spectrum is 12 A. Indeed it was included in the 
list of ISIS service targets because of the suspicion that this 
might be the case: when it is plotted on the colour-colour 
plane with only neighbouring stars, within a 10 x 10 arcmin 2 
box, it sits clearly separated in (r — Ha) just above the 
local main stellar locus (see Fig. 12). This is an example of 
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Figure 12. The (r' — Ha, r' — i') plane in the locality of star 
B. The plotted objects are selected with 13 < r' < 19 and are 
located in a 10 X 10 arcmin 2 sky area centred on RA 05 33 40 Dec 
+25 20 00 (2000). In the absence of its weak Ha emission, star B 
would have fallen 0.13 lower in (r 1 — Ha), at a position typical of 
non-emission A0V stars. 

the greater coherence of mean conditions (reddening, nature 
of population) within a smaller sky area leading to greater 
success in identifying an 'unusual' object. 



6 SPECTROSCOPIC TRAWLING FOR 
EMISSION-LINE AND OTHER RARE 
OBJECTS - FIELDS IN CEPHEUS 

We now present some results from early spectroscopic follow- 
up of IPHAS in order to give a concrete example of the yields 
of different object types from IPHAS data. In this case, the 
choice of sky area has been dictated mainly by observational 
convenience, rather than by data quality considerations. The 
results presented here rest on more typical IPHAS photom- 
etry. 

6.1 The MMT/HectoSpec observations 

In June 2004 we obtained spectra with the Mount Hopkins 
6.5-metre MMT in F/5 configuration using the recently com- 
missioned HectoSpec facility, a multi-object spectrograph 
fitted with 300 fibres that can be deployed across a field, 
1 degree in diameter (Fabricant et al 2004). The fibre posi- 
tioner is mounted at Cassegrain. The 270 groove/mm grat- 
ing used delivers broad wavelength coverage (4488 - 8664 A) 
at 6.2 A resolution. Over two nights, six different fields were 
observed using two fibre configurations per pointing. The 
target stars selected for this programme fell mainly in the 
magnitude range 17 < r' < 20. The total on-source exposure 
times were 1200 sees. Spectra were extracted by the instru- 
ment pipeline that includes CCD bias and gain corrections, 
flat-fielding using domeflats as well as a sensitivity correc- 
tion for the individual fibres using twilight flats. Individual 
fibre spectra were then extracted and wavelength calibrated 
using FeNeAr-lamp exposures. Finally, a mean sky spectrum 
derived from the sky fibres was subtracted. Due to spatial 
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Figure 13. Point-source colour-colour data derived from IPHAS 
fields 7012/7012o in Cepheus. This plot illustrates the quality 
of data around (I = 105.6°, b = +4.0°) from which we selected 
targets for multi-fibre spectroscopy with MMT/HectoSpec. The 
magnitude range from which HcctoSpec targets were selected was 
(17 < r' < 20) (blue data points in both panels - data outside this 
magnitude range are plotted as smaller black points). The black 
lines in both panels are the simulated unreddened main sequence 
track and early-A reddening line, shifted downward in (r' — Ha) 
by 0.16 magnitudes. Here they serve as fiducial lines to show how 
the errors grow at fainter magnitudes. 



variations in the sky background across the field of view, sky 
subtraction using sky fibres is not always perfect and care 
must be taken with spectra displaying weak, unresolved Ha 
emission components. Sky subtraction will be improved in 
the future using offset sky exposures that sample the sky for 
each fibre close to the target position, in conjunction with 
an optimal scaling correction using the strongest sky lines. 

6.2 Selecting targets for the Cepheus field 

We report the results from one of two pointings in 
the constellation of Cepheus. Centred on RA 22 17 
00, Dec +61 33 37 (2000) {I = 105.6°, b = 4.0°), 
this position was picked because it contains a strip 
included in the Spitzer Galactic First Look Survey 
(http://ssc.spitzer.caltech.edu/fls/galac/). Despite its rela- 
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Figure 14. The point-source colour-colour diagram used for 
target selection in the l°-diameter field in Cepheus observed 
with MMT/HectoSpec. Data for objects in the magnitude ranges 
17 < r' < 20 and 20 < r' < 20.5 are shown in blue and black 
respectively. Because of IPHAS field merging and the accompa- 
nying ad hoc photometric corrections, the main stellar locus is 
more diffuse here than in individual field extractions (cf Fig. 13). 
The simulated unreddened main sequence track and the early-A 
reddening line (drawn as solid lines) are accordingly less useful in 
decoding the data. 



tively high Galactic latitude, this area of sky presents sig- 
nificant and locally-variable interstellar extinction (ranging 
from Av ~ 4 up to ~7 magnitudes). This shows up in 
the IPHAS colour data for the region as somewhat broad- 
ened main stellar loci. An example is shown as Fig. 13, 
where colour data extracted from the IPHAS field pair, 
7012/7012o, are plotted. At magnitudes brighter than r' ~ 
19 (top panel), the colour uncertainties are less than ~0.05 
and are less significant than environmental factors in the 
smearing of the main stellar locus - this reverses at fainter 
magnitudes (lower panel) where the typical errors are ~ 0.1. 

The HectoSpec l°-diameter field spans then IPHAS 
field positions, not including offsets. The largest contribu- 
tions, however, are from IPHAS fields 6985, 6993, 7012 and 
7019 (see Table 1). To speed up the compilation of the target 
lists, we chose to merge the data from all the relevant point- 
ings first, before proceeding to target selection. In merging 
the data, corrections for photometric shifts between different 
WFC exposures had to be applied. These were calculated us- 
ing the mean magnitude offsets for sources located in field 
overlaps. Inevitably, this merging blurred the main stellar 
locus in the colour plane some more - compare the lower 
panel in Fig. 13 with the plot for the full HectoSpec field in 
Fig. 14. For the future, we are re-ordering the algorithm, in 
order to give more emphasis to selection at the individual 
field level (cf the discussion of 'star B' at the end of Section 
5). Indeed a further tactic that can be applied in order to 
minimise the spread of the main stellar locus is to select 
from within a number of narrowly-set r' magnitude ranges. 
On this occasion, after merging the catalogues for the rele- 
vant individual fields together, the main target selection was 
performed within the full magnitude range to be observed, 
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17 < r' < 20. Finally a few promising emission line star 
candidates down to r' = 20.5 were added by hand. 

The goal of this first round of HectoSpec observations 
was to explore the complete IPHAS colour-colour plane 
while trying to give high priority to objects that are outliers 
or near the edge of the general distribution of objects. This 
naturally includes all emission line star candidates, which 
lie above the main stellar locus. 

To achieve this sampling, the following selection algo- 
rithm was applied. The colour-colour plane was split into 
bins of 0.1 magnitudes in width and height. Targets were 
then selected based on the number of objects in each bin: 
from 1 to 3 objects in a bin, all were selected and given the 
highest fibre allocation priority; between 4 to 9 objects in a 
bin, a random selection of 75% of the bin members was cho- 
sen with a slightly lower fibre allocation priority; between 10 
and 50 objects, a random fraction falling linearly from 50% 
to 10% was selected and given a low fibre selection priority 
based on the number of objects in the bin; >50 objects in a 
box, 10% of the objects were selected randomly and given a 
low fibre selection priority based on the number of objects 
in the bin. The maximum number of objects that could be 
selected in a box was capped at 10. Finally, if a bin had 4 
neighbouring bins in the cardinal directions with 10 or more 
objects, then only one object was selected for spectroscopy. 
This rule was introduced to further reduce the number of 
objects picked for spectroscopic follow up that lie in the 
densest - and probably most uninteresting - parts of the 
stellar distribution. 

If the selection algorithm gave less than the requested 
number of objects (~ 700 for 2 HectoSpec configurations) 
then 1 object per bin was added, starting with the least 
populated bins until the required number of objects was 
reached. If more than the requested number of objects was 
chosen, then 1 object was removed from each bin with more 
than 1 object selected - starting with the most populated 
bin until the requested number of objects was reached. By 
this means we were able to assign fibres to between 450 and 
540 targets, depending on the field observed. Each set was 
then split into two configurations. 

The outcome of this process and the HectoSpec obser- 
vation of the Cepheus field was a collection of 496 stellar 
spectra. 

6.3 Results of the MMT/HectoSpec spectroscopy 
in Cepheus 

First and foremost we find that essentially every target star 
located in the colour-colour plane clearly above the main 
stellar locus is confirmed as an emission line star. Altogether 
29 objects are confirmed as having Ha in emission, with 6 
or 7 lying on or just below the upper bound of the highly 
populated region (green encircled points in Fig. 15). The 
one object not encircled in green, at (r' — i') ~ 1.1 and 
(r — Ha) ~ 1.5, is probably an emission line object also, 
but remains ambiguous because it is very faint and its spec- 
trum is correspondingly noisy. There is a still larger group, 
numbering 47, of probable dMe stars. They are 'probable' 
because the sky subtraction may have left a false residue of 
Ha emission. The range of Ha emission equivalent widths 
in this group is from a few up to 10-20A. Quite plausibly, 
most of these objects lie mixed in with non-emission dwarf 



M stars (of which there are 90 or more). The one 'probable 
dMe' stars (at (r' - i') ~ 2.50, (r' - Ha) ~ 1.3) above the 
main locus can be viewed as dMe with the greatest confi- 
dence because time-variable Ha emission is characteristic of 
dMe stars - presumably at the time of the IPHAS imaging, 
its Ha emission was brighter than 9 months later, at the 
time of the spectroscopy. 

The combination of moderate spectral resolution and 
short exposure time has meant that many of the more rou- 
tine objects, without Ha in emission or in marked absorp- 
tion, are more challenging to sort into spectral classes. This 
large group of 214 objects will be dominated by late- A to 
mid K stars, but will also include some non-emission OB 
stars. The spectra of a further 89 stars are so faint and 
noisy that no comment can be made about them. The stand- 
out objects towards the lower boundary of the main stellar 
locus in Fig. 15 are the stars with Ha strongly in absorp- 
tion. There are 15 of these. Two of them are well-separated 
from the main locus at much lower (r' — i') and also much 
lower (r' — Ha) - they are both white dwarfs. Similarly 
placed objects in other IPHAS fields for which we have 
MMT/hectospec spectra have turned out to be white dwarfs 
too. The remaining thirteen stars with strong Ha absorption 
are early A stars. Tighter classification at this time is not 
feasible. 

We now present a selection of 8 objects and their spec- 
tra for more detailed discussion. These are identified in the 
colour-colour plane shown as Fig. 15. The data on them, 
given in Table 6, includes estimates of their r' magnitudes 
and (r' — i'), (r' — Ha) colours derived from the highest 
quality IPHAS exposures currently available. Note that the 
colours are, typically, different from those plotted in Figs. 15 
- this is due to the colour shifts applied in combining IPHAS 
fields before MMT/HectoSpec target selection. Since they 
assist in assigning broad object class, we also include in the 
table 2MASS (J - H), (H - K) colours and K magnitudes. 

Stars 1-3 (Fig. 16) are most likely to be young stellar 
objects (YSOs) of Herbig or T Tau type. This object class 
assignment is easiest for star 1 since the veiling is not so 
extreme as to hide the underlying M-star spectrum. Indeed, 
a comparison between stars 1 and 4 in Fig. 16 suggests that 
these objects' M spectral sub-types are likely to be very sim- 
ilar. The sky around star 1 has been imaged in all four IRAC 
bands by the Spitzer Space Telescope First Look Survey. We 
downloaded the calibrated images from the Spitzer Science 
Archive and carried out point source extractions. Star 1 was 
detected at 3.6, 4.5 and 5.8 /jxa, with fluxes corresponding 
to magnitudes of 12.47±0.11, 12.12±0.12 and 11.93±0.14, 
respectively. Using the observed IPHAS, 2MASS and IRAC 
fluxes, we find that the SED of Star 1 from the r'-band to 
5.8 fim happens to fit quite well to a 2000 K blackbody (al- 
though the SED must include components due to a reddened 
stellar photosphere, an accretion disk and warm dust). 

Star 2 has the highest contrast and richest emission line 
spectrum of the three stars - in this object, even the Na I 
D lines are thrown into emission and some forbidden line 
emission is present. Star 3 is intermediate between 1 and 
2, both in terms of the contrast of its emission spectrum, 
and that it is just possible to pick up late-type photospheric 
absorption against the continuum (e.g. the blend at 6495 A 
strong in G/early-K stars). The NIR colours of all 3 objects 
imply modest NIR continuum excesses that are not out of 
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Figure 15. The colour-colour diagram of the Cepheus targets observed with MMT/HectoSpec. Symbols used: each small black circle 
locates the IPHAS colours of an object for which we have a spectrum; larger green circles pick out confirmed emission line stars; red 
circles pick out weaker-emission, possible dMe stars; larger black circles indicate non-emission line late-type stars with strong molecular 
bands in their spectra; blue circles pick out stars with prominent Ha and H/3 absorption features; black crosses mark those objects for 
which the spectrum is too noisy to classify. Small black circles left unadorned will mainly represent late-A to early K stars. The synthetic 
tracks, drawn as solid lines are as in Figs. 13 and 14. Spectra of the objects 1—8 are shown in Figs. 16 amd 18 and parameters on them 
are given in Table 6. 



Table 6. Positions, magnitudes and colours for the 8 stars, observed using MMT/hectospec, whose spectra are plotted in Figs. 16 and 
18. The r' magnitudes and (r' — Ha), (r' — i 1 ) colours have been taken from the catalogues for the best pair of IPHAS exposures. At this 
time, the magnitudes are likely to be correct to within ±0.1, while the Ha EWs are reliable to within ±5 A. The EW sign convention is 
reversed in that a positive value implies net emission. 





IPHAS name/position 




IPHAS photometry 


object 


2MASS colours/maj 


mitudes 


Ha EW 




J[RA(2000)+Dec(2000)] 


r' 


r'-i' 


r' - Ha 


type 


(J-H) 


(H-K) 


K 


(A) 


1 


J221734.39+611409.2 


19.5 


1.55±0.06 


1.72±0.06 


T Tau star 


1.22±0.07 


0.67±0.05 


13.39±0.03 


190 


2 


J221740.30+614702.9 


17.7 


1.44±0.01 


1.42±0.01 


B[e]/YSO 


1.33±0.04 


0.98±0.04 


10.52±0.02 


200 


3 


J221411. 60+612606.7 


17.5 


1.15±0.01 


0.99±0.01 


Be/T Tau 


1.38±0.08 


0.92±0.05 


12.12±0.03 


80 


1 


J221427.30+612943.6 


19.5 


2.01±0.08 


1.08±0.08 


~M3Ve 


1.09±0.04 


0.56±0.04 


13.15±0.02 


70 


5 


J221534.47+615725.4 


17.2 


0.06±0.01 


-0.23±0.01 


white dwarf 








-34 


6 


J221345.86+614418.6 


17.5 


0.53±0.01 


0.00±0.01 


early A star 


0.3±0.2 


-0.5 


~15.3 


-16 


7 


J221822.09+614803.8 


18.7 


3.39±0.03 


0.78±0.04 


mid-M giant 


1.66±0.04 


0.62±0.04 


8.68±0.02 




8 


J221619.89+612621.1 


19.5 


2.97±0.06 


0.71±0.07 


carbon star 


2.61±0.03 


1.76±0.03 


6.41±0.02 
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place for Herbig or T Tau stars, with star 2 showing the 
most marked excess. 

It seems likely that these stars 1-3 are associated with 
LDN 1188 (£ = 105.7, b = +4.2), a dark cloud less than 2° 
away from the well-known Sh 2-140 region of star formation. 
This proximity suggests a similar distance to both regions, 
which has been given as 910 pc for Sh 2-140 by Crampton 
& Fisher (1974). However both these nebulae lie on the pe- 
riphery of Cep OB2, to which the distance appears to be 
rather less (~600 pc, de Zeeuw et al 1999). In a study of 
LDN 1188, Abraham et al (1995) reported the discovery of 
a number of emission line stars in objective prism data ob- 
tained at Konkoly Observatory at brighter magnitudes than 
our HectoSpec selection. RNO 140 and RNO 141 (Cohen 
1980) are also in this neighbourhood. 

Star 4 can be presumed to be in the foreground with 
respect to LDN 1188 and stars 1-3. Hence its JHK colours, 
combined with its spectral type can be used to estimate 
a minimum interstellar extinction towards LDN 1188. This 
works out at E(B - V) ~ 1.3 or A v ~ 4 (for R = 3.1, 
and using data from Bessell & Brett 1988). We can now 
see if this marries up with the implications of the IPHAS 
colours of stars 1 - 3 by comparing their catalogue values 
with reddened synthetic estimates. This is accomplished via 
Fig. 17 in which the IPHAS colours for stars 1 to 3 (Table 6) 
are compared with synthetic tracks (Table 4) . 

If star 2, with its generally high contrast emission line 
spectrum is a Herbig or T Tau star with an accretion domi- 
nated SED, its reddening would correspond to E(B— V) ~ 2: 
for, in Fig. 17, it lies just to the right of the track for 
Fx oc A~ 2 ' 3 and E(B — V) — 2. A similar, or somewhat 
lower, reddening would appear plausible for star 1 in that its 
optical SED should be somewhat redder, intrinsically, than 
that of star 2. The intrinsic optical SED of star 3 should 
be intermediate between stars 1 and 2 (given the marginal 
detection of late- type photospheric absorption), and yet it 
is observed to be 'bluer' than either. The highest likely red- 
dening of star 3 is E(B — V) ~ 1.5: this reddening would 
apply in the limiting case of an accretion-dominated SED, 
where the contribution of the G/K star is small. 

On the basis of its (r' — i') colour and spectral type, and 
after correction for its Ha emission, star 4 would be assigned 
E(B — V) ~ 1. This is a somewhat lower estimate than the 
estimate based on 2MASS NIR colours (E(B - V) ~ 1.3) 
but not so large a discrepancy that either the NIR or optical 
photometrically calibration must be called into question. In 
conclusion, we find that the order of increasing reddening 
appears to be: star 4 in the foreground, star 3, and then 
star 1 and star 2, spanning the range 1 < E(B — V) < 2. 
The Schlegel et al (1998) Galactic reddening maps indicate 
maximal reddenings of E(B — V) ~ 2.5 for this part of the 
Plane. We have rough consistency and a first indication of 
patchy reddening toward the young objects in the vicinity 
of LDN 1188 . If half the r' flux of star 1 is attributed to 
an ~M3 stellar photosphere, and Av ~ 6, one may deduce 
a stellar radius of around 4 times the M3V main sequence 
radius, for a distance of ~ 600 pc. 

The IPHAS colours for stars 2 and 3, picked out in 
Fig. 17, are broadly consistent with their observed Ha EWs 
(200 A and 80 A), in that their colours 'predict' EWs of 
~200 A and ~100 A, respectively. Only star 1 is discrepant 
in that its observed EW (~190 A) is distinctly low com- 
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Figure 16. The MMT/HcctoSpcc counts spectra of Cepheus tar- 
gets 1 to 4, as marked in Fig. 15 and listed in table 6. These data 
and those shown in Fig. 18 are uncorrected for telluric absorption. 



pared with the implications of its IPHAS colours (suggest- 
ing an enormous EW of about 250 A). Given that Ha EW 
is well known to be a time-variable quantity in most classes 
of emission line object, consistency for 2 out of 3 objects is 
acceptable. 

The spectra of 4 non emission line objects, star 5-8, 
are shown in Fig. 18. These draw attention to what can be 
called Ha deficit positions in the IPHAS colour-colour plane 
(cf Fig. 15). At the blue end, objects 5 and 6 are examples 
of extreme and strong Ha absorption objects: respectively a 
white dwarf and an early A star. At the red end, objects 7 
and 8 are respectively a normal, somewhat reddened mid-M 
giant, while object 8 has the distinctive CN band structure of 
a carbon star in its spectrum. Object 7 is typical of the stars 
populating the red end of the giant strip in the IPHAS plane. 
The NIR colours and rough spectral type (~M4III) suggest 
a reddening corresponding to E(B — V) of about 2. Carbon 
stars like object 8 will usually fall below the red giant strip 
(lacking the TiO bands that, in M-type spectra, lead to the 
seeming flux maximum in the Ha region) . In terms of its NIR 
colours, star 8 is a more extreme object than the suspected 
carbon star mentioned at the end of section 5 (although 
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Figure 17. A diagnostic colour-colour diagram for the Cepheus 
region stars 1 to 3 in Table 6. The blue, light blue, (right hand) 
green and red tracks, from Table 4 and Fig. 6, are reddened to 
E(B — V) = 2. They respectively represent the colour-colour re- 
sponse to increasingly strong (narrow) Ho line emission of ex- 
treme O stars, ~A0 stars, optically-thick accretion disks and ~G2 
stars. For contrast, the track for optically thick disk accretion at 
E(B — V) = 1 (left hand green line) is also shown. The hor- 
izontal dashed lines are lines of constant Ha equivalent width. 
The numbers 1-3 mark the positions of stars 1-3: note that their 
(r' — Ha) colours have been incremented by 0.16 (cf Fig. 13) to 
place them in the simulation colour domain. The required correc- 
tion in (r' — i') will be significantly smaller. 



reddening gives both objects more extreme (J — H , H — K) 
colours than seen in the local Galactic carbon star sample 
of Claussen et al, 1987). Indeed Object 8 appears to be a 
reddened (A v ~ 5) version of 2MASSI J0326599+ 143957, 
described by Liebert et al (2000) as a luminous, very cool, 
late N type carbon star. 

In summary, our early MMT/Hectospec observations 
of this field in Cepheus sample a diverse range of objects. 
Sources lying clearly above the stellar locus in the (r — 
Ha, r' — i') colour-colour plane have indeed been confirmed 
as true emission line objects with Ha EWs ranging from 
a few to 200 A. Many of them are likely to be Herbig or 
T Tau stars. In addition several Ha deficit sources have been 
identified. Multi-object spectroscopic follow-up will remain 
a key part of our efforts to mine the IPHAS database. 



7 IPHAS OPPORTUNITIES FOR Ha IMAGING 

The power of the survey for detecting stellar Ha emission 
has already been described. A second target for the survey 
are the spatially-resolved emission-line nebulae. These neb- 
ulae are - like emission-line stars - associated with early 
and late phases of stellar evolution, and indicate ionization 
of circumstellar gas, in the form of H n regions, planetary 
nebulae and supernova remnants. Already, in the southern 
hemisphere, the UKST SHS has shown the remarkable in- 
completeness of existing catalogues by doubling the number 
of known planetary nebulae (Parker et al 2003). 

The Galactic Plane shows ubiquitous diffuse Ha emis- 
sion as well as reflection nebulae, which need to be distin- 
guished from circumstellar ionized nebulae. Since reflection 
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Figure 18. The MMT/HectoSpec counts spectra of Cepheus tar- 
gets 5 to 8, as marked in Fig. 15 and listed in Table 6. 



nebulae are the product of continuum scattering, they can 
be removed by comparison of Ha and r' images. Diffuse 
Ha emission occurs on large size scales (10 arcminutes to 
degrees) and lacks the usual symmetry of circumstellar neb- 
ulae. 

Imaging of extended nebulae requires well behaved 
background on the CCDs making up the WFC. Because of 
this, the preferred image properties for nebular studies are 
very different from those for point source extraction: high 
background is acceptable for the latter, while poor seeing is 
acceptable for the former. Small nebulae covered well within 
a single CCD do not require special reduction. Each object 
will usually have been covered in at least two pointings by 
the time the survey is complete, giving improved S/N. For 
larger objects, mosaics need to be made. Subtracting an r' 
frame removes stars: for small areas this can be done using 
PSF matching techniques, but this is very computationally 
intensive - such that for larger fields, a direct subtraction is 
used which typically leaves larger residuals. 

The limitations of the technique are largely due to back- 
ground variations. The observations are generally taken in 
grey and bright time. Different fields will therefore present 
with very different background sky levels. The sky subtracts 
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Figure 19. The continuum-subtracted Ha image mosaic of the extremely large angular size supernova remnant, S 147. The frame size 
is roughly 5 degrees across, and 4 degrees down. The extent of the SNR is given by Van den Bergh et al (1973) as 200 x 180 arcmin 2 . N 
is up and E to the left. The H II region, Sh 2-242 is picked out near the eastern edge of the image by the white box (see Fig. 20). 



fairly well in an Ha — r' image (unless the background is 
variable under non-photometric conditions), but if there is 
smooth extended Ha emission over large angular scales, its 
contribution is currently not separately determined from the 
sky contribution. Internal reflections are seen in some im- 
ages, from bright stars. In some locations, a bright star just 
outside the field of view gives a flare-like feature on the edge 
of a nearby frame. 

A photometric calibration is determined for point 
sources, but not for extended emission. To correctly calibrate 
emission nebulae, assuming the continuum background is 
fully subtracted, the filter response curve needs to be pre- 
cisely known and needs to be stable over the likely Doppler 
wavelength shifts (Ruffle et al 2004) . [N n] will also intrude 
into the flux. Currently, calibration is best performed using 
known planetary nebulae located in the imaged area. Note 
that stellar Ha sources cannot be used for calibration as 
their line flux tends to be time variable. 



7.1 A supernova remnant 

As an example of the possiblities regarding extended 
nebulae, we present an image of the supernova rem- 
nant S 147 (Shajn 147, or Simeis 147, in full - not 



Sh 2-147, with which it is confused in SIMBAD). This 
is a near-perfect remnant of an approximately spheri- 
cal shape, showing a typical filamentary structure. It is 
positioned just overlapping the anti-galactic centre, its 
own centre being at £ ~ 180.1°, b ~ —1.6°. But due 
to its large extent, spanning several degrees, only pho- 
tographic images have been published so far (Van den 
Bergh, Marscher & Terzian, 1973). See also the 24th March 
2005 'Astronomy Picture of the Day' due to R. Gendler 
(http:/ /antwrp. gsfc.nasa.gov/apod/ap050324. html). 

An IPHAS image of S 147 was produced by combining 
approximately 250 pointings. The pipeline mosaicing proce- 
dure was found to be inadequate for combining many fields 
taken under widely varying conditions. We therefore used 
the pipeline only to produce reduced images of the indi- 
vidual CCDs for each pointing. Then, for each image, we 
subtracted the r' image from the Ha data and smoothed to 
a pixel size of 5 arcsec (a binning factor of 15). The back- 
ground per image was approximately nulled by subtraction 
of the median: in fields with bright and extended emission 
this required manual selection of areas for background def- 
inition - otherwise the median over all four CCDs of one 
pointing was used. One corner of CCD no. 3 is affected by 
scattered light in conditions of bright Moon light (amplified 
if cirrus was present): this corner was always blanked out. 



Ha survey of the northern Galactic Plane (IPHAS) 23 




Figure 20. The Hn region, Sh 2-242, picked out on the eastern 
edge of Fig. 19, shown at full resolution. This is a false colour 
composite of the images in the three filters. The imaged area is 
roughly 10 x 10 arcmin 2 . N is up and E to the left. 



All CCDs tend to show a gradient along the long axis of 
about 1 ADU in the subtracted image: this could be due to 
a charge-transfer efficiency limitation, but its exact cause is 
not known. It does not subtract automatically because the 
r' image is scaled first - accordingly, as a final step, this 
remaining gradient was subtracted from all images. 

The resulting images were combined in a single mosaic 
using the Virtual Observatory software package Montage. 
This involved regridding each frame, and determining back- 
ground corrections by comparing areas in common between 
different images. We finally produced an image covering 25 
square degrees. We note that extended emission on scales 
of a degree or more may not be well represented, as it can 
be affected by the background subtraction procedure. How- 
ever, the filamentary structures are very well recovered. The 
background fitting was found to be insufficiently constrained 
at the outer edges of the imaged area, leaving some nega- 
tive areas. To deal with this, a linear gradient was fitted 
to the background in empty regions of the full field, and 
subtracted. 

The image of S 147 is shown as Fig. 19. A bigger ver- 
sion is to be found in the INC Newletter article describing 
IPHAS (Drew et al 2005). The improvement is very marked 
in comparison with the photographic image shown by Van 
den Bergh et al (1973): the huge increase in dynamic range 
brings with it a subtlety of detail missing from the old im- 
agery. Furthermore, these data allow smaller areas to be 
imaged to much higher resolution. But even at this 5 arcsec 
rebinning the structure is clear. We note that a blow-out is 
obvious on two sides of the remnant: left and right (to the 
E and W) . The full extent of S 147 is essentially as reported 
by Van den Bergh et al, at a mean diameter of just over 3° 
- several times that of the Moon. 

At the eastern edge of the image mosaic, and south 
of centre, lies another prominent but much more compact 
nebulosity, Sh 2-242. This very bright smudge in Fig. 19 has 
a mean Ha diameter of 8.0 arcmin. As a contrast to the 



very large scale structure of S 147, we show as Fig. 20, a full 
resolution image of this H II region. 



8 SUMMARISING DISCUSSION 

The main aim of this paper has been to introduce IPHAS, 
the INT Photometric Ha Survey of the Northern Galactic 
Plane. By the beginning of 2005, 55% of the imaging obser- 
vations had been obtained. It is expected that the survey 
will reach completion in 2006. 

Given that previous surveying of the northern Galactic 
Plane for emission line objects rarely reached deeper than 
V ~ 13 (see KW99), while the sensitivity limit of IPHAS 
is r' ~ 20, there is no doubt that a huge domain is being 
opened up for exploration for the first time. It is difficult 
to predict the numbers of emission line objects that will be 
discovered, even now with the survey in progress, because 
the distribution of such objects along the Galactic Plane is 
extremely uneven (see KW99). The small number of IPHAS 
fields discussed in this paper fit in with this impression: 
no emission line objects were evident in the colour-colour 
plane for the Taurus field (2540/2540o in section 5 - but 
a very modest example was found on closer examination); 
a handful were evident in the Aquila fields (section 4), and 
upwards of 20, all fainter than r' — 17, were picked up in 
the MMT/HectoSpec pointing toward LDN 1188 in Cepheus 
(section 6). Crudely averaging the experience to date, it is 
likely IPHAS will uncover around 10 emission line objects 
per square degree in the range 13 < r' < 20 (roughly 2 to 3 
per IPHAS field), and hence no less than ~20,000 altogether. 

Through its use of both narrow-band Ha and broad- 
band r' and i' CCD photometry, IPHAS has the capability 
to pick out Ha deficit objects - a possibility typically be- 
yond objective prism spectroscopy, the traditional tool of 
emission line star hunting. For example, unreddened white 
dwarfs are easily spotted as blue objects with small or neg- 
ative (r' — Ha), separated from the main stellar locus in 
the IPHAS (r' — Ha, r' — i') plane. It is a reasonable guess 
that over a thousand will be discovered in IPHAS data, just 
as have been discovered in a comparable sky area around 
the north Galactic cap via the Sloan survey (Kleinman et 
al 2004). Very red Ha deficit objects can be either brown 
dwarfs or carbon stars - classes of star that both lack the 
TiO band absorption responsible for raising (r 1 — Ha) in nor- 
mal M stars - or they can be very reddened, rare examples 
of late-type supergiants. Two carbon stars (one probable, 
the other confirmed) have been reported here. 

A further role for IPHAS is that it can help trace the 
way in which the stellar populations making up the Galactic 
Plane vary across the northern sky. It has been shown here 
that not all sightlines look the same in IPHAS colours: of 
particular note is the sharp contrast between the red-giant 
deficient Taurus field (£ — 181.7°, section 5), in the Galactic 
Anticentre region, versus the Aquila fields (£ ~ 33°) with 
their prominent giant-star populations, sampling the inner 
Galaxy. Before now, rather little has been known about the 
far reaches of the Milky Way outside the Solar Circle. For 
example, the recent study of Galactic spiral structure by 
Russeil (2003) reaches only to ~ 6 kpc, outward from the 
Sun. The sensitivity of IPHAS is more than adequate for ex- 
tending our knowledge much further: even relatively humble 
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A0-A3 V stars - easily picked out around what has been 
dubbed the early-A reddening line in the IPHAS colour- 
colour plane - are potentially accessible out to distances of 
~20 kpc in the direction of the lightly reddened Galactic 
Anticentre. 

We have laid out the character of the (r' — Ha, r' — i') 
colour-colour plane that is unique to this survey, and have 
established a grid of simulated colours that will be of use 
in the analysis of IPHAS observations. In doing this, we 
have inevitably identified items of work for the future. For 
example, a job remains to be done to achieve as good a 
quantitative match between simulations and observations of 
M stars as is achieved for earlier spectral types. A related 
problem has been noted by those developing Sloan Digital 
Sky Survey photometric calibrations (Smith et al 2002) . The 
reward for solving this problem will be the chance to fully 
realise the potential of IPHAS as an unparalleled resource 
for the statistical analysis of M-dwarf activity. Another task 
is to gather a a range of early A star spectra that can be 
used to achieve better definition of the early-A reddening 
line and its dependences. As specified here, it is unlikely to 
be more than ~0.02 magnitudes from its correct registration 
in {r' -Ha). 

Bigger issues to be dealt with are that the photomet- 
ric calibration of IPHAS data will need to be made uniform 
across the survey, and that a proper definition of Ha mag- 
nitude zero point (referred to Vega) is needed that can be 
applied within the pipeline processing. Until these calibra- 
tion requirements have been met, it will remain necessary to 
derive colour offsets between observed and simulated data 
independently for every IPHAS field. Typically, offsets in 
(r 1 — i') will be small - but for (r' — Ha), the differences be- 
tween catalogue and simulated colours may vary from ~0.1 
up to ~0.2 magnitudes. In most cases it is easy to gauge 
the required offset by matching the synthesised unreddened 
main sequence track and early-A reddening line to the up- 
per and lower boundaries of the main stellar locus in the 
colour-colour plane. 

Simulations have also been performed that show how 
the equivalent width threshold for the detection of Ha emis- 
sion will change with observed [r' — i') colour (section 3.2 
and Fig. 6). The photometric accuracy of IPHAS is such 
that (r' — Ha) differences of 0.05-0.1 are significant down 
to r' ~ 20. Re-expressed in terms of an Ha equivalent width, 
and for magnitudes brighter than ~19, this corresponds to 
a threshold emission EW of roug hly 5 A. 

It has been shown that the typical morphology of the 
main stellar locus in the colour-colour plane permits the 
selection of candidate emission line stars presenting with 
threshold Ha emission at low (r' — i') only. In practise the 
IPHAS bright magnitude limit (r' ~ 13) has the effect of 
all but eliminating normal stars with colours bluer than 
(r — i') ~ 0.5 from the point source catalogue (note Figs. 3, 
9 and 14) - with the consequence that lightly-reddened sub- 
luminous accreting objects, with or without Ha emission, 
will usually lie comfortably outside the main stellar locus. 
IPHAS can therefore be used straightforwardly to identify 
all such objects (in addition to many non-interacting white 
dwarfs, found as 'deficit' objects). 

The threshold for Ha emission high-confidence detec- 
tion rises from ~ 10 A equivalent width at (r' — i') ~ 1 up 
to ~ 50 A at (r' — i') ~ 2.5 - beyond this the colour-colour 



plane typically becomes sparsely populated again. This has 
important implications for the detectability of T Tau and 
other young emission line objects: based on equivalent width 
data collected by Reipurth, Pedrosa & Lago (1996, their Ta- 
ble 1 and Fig. 10), one-third to a half of such objects would 
be immediately identifiable as emission line objects in the 
IPHAS database at E(B - V) ~ 3. At lower reddenings 
the fraction would be higher. This suggests that the roles 
for IPHAS with respect to young stellar populations are (i) 
finding extreme examples anywhere (here we have presented 
two with Ha emission EWs of ~ 200 A), and (ii) picking 
out new associations through the detection of its most ac- 
tive members. Once a new association is identified, a more 
detailed exploration of catalogued IPHAS sources, exploit- 
ing the tricks of e.g. apparent magnitude binning, or small 
area searches, would be likely to uncover further candidate 
emission line sources. 

Evolved high mass stars, with Ha emission, that find 
their way into the IPHAS database will do so because they 
are very distant and significantly reddened (they are other- 
wise too bright). The more extreme, least well-understood, 
and therefore more interesting groups (Wolf-Rayet stars, lu- 
minous blue variables, B[e] stars, yellow hypergiants) usually 
present with very high EW Ha emission (~ 100 A and more) 
and so will not be missed. The most extreme point-source 
emission line objects of all, compact nebulae, are most likely 
to appear to be relatively blue in (r' — i') or even evade de- 
tection in the i' band, and will have the most extreme values 
of (r' — Ha) possible, i.e. not more than 3 to 3.1, in practice. 

This connects naturally to a topic only touched on in 
this paper - the exploitation of IPHAS in the study of 
spatially- resolved nebulae (deferred to a later paper). A 
search for PNe is in progress, in which several tens of can- 
didates have been identified and a few have been studied 
spectroscopically. A paper on the study of an intriguing 
quadrupolar nebula located well-outside the Solar Circle is 
in preparation (Mampaso et al) . Here we have simply drawn 
attention to the flexibility IPHAS presents for the investiga- 
tion of a wide range of spatial scales, ranging from arcsec- 
onds to several degrees. 

Further power to diagnose either particular object types 
or complete stellar populations will come from pooling 
IPHAS r' , Ha and i' data with data from surveys in differ- 
ent wavebands. In the future of an astronomy conducted via 
virtual observatories, a survey as comprehensive as IPHAS 
- with its particular exploitation of narrowband Ha data - 
will be a major resource. Currently there is an obvious syn- 
ergy with the all-sky NIR survey 2MASS, although this is 
limited to the reddened parts of the northern Galactic Plane 
since 2MASS reaches only to K ~ 15. The gap this leaves 
should soon be plugged by the UKIDSS Galactic Plane Sur- 
vey, reaching to K — 19 (see http://www.ukidss.org/). Be- 
yond IPHAS, comprehensive optical surveying of the north- 
ern Galaxy using linear detectors is still lacking. However it 
is interesting to note the recent release by SDSS of u'g'r'i'z' 
data on a number of low Galactic latitude fields (Finkbeiner 
et al 2004). 

We finish with a comment on the plans for making 
IPHAS data available to the community. At the present 
time there is open access to the reduced images held at 
CASU, from a year after the images have been processed. 
Accordingly the images obtained in 2003 can already be 
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accessed by anyone able to reach the CASU web inter- 
face (http://apm2.ast.cam.ac.uk/cgi-bin/wfs/dqc.cgi). Ac- 
cess immediately after processing is available to those work- 
ing in the ING partner countries: the United Kingdom; 
Spain and The Netherlands. The point source catalogues 
are to be released in two stages: in the first half of 2006 we 
aim to release as many as are available, calibrated only at 
the individual exposure level (as they have been described 
here); we intend to follow this up with a second release, after 
the survey is complete, and when a uniform calibration has 
been established across all fields. We anticipate that the fi- 
nal catalogue will contain photometry on around 80 million 
point sources, capturing data on roughly 1 for every 1000 
stars estimated as existing in the northern Milky Way. 
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